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1. Introduction 
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Gamma-ray bursts (GRBs) and supernovae (SNe) bring new perspectives to the study 
rf-"] ^ of neutron stars and white dwarfs, as well as opening new branches of theoretical physics 

y—^ I and astrophysics. 

^' 

Q , I dedicate this talk to John Archibald Wheeler recalling some crucial moments in our 

M ' collaboration and focusing on two of the most energetic and transient phenomena in 

^ ■ the universe: supernovae (SNe) and gamma-ray bursts (GRBs). From the knowledge 

being gained daily on these astrophysical systems a new physical and astrophysical 
understanding is emerging. The physics of neutron stars and black holes is rapidly 
^ , evolving, new more complex astrophysical scenarios previously unimaginable are 

^SJ I being developed and new domains of fundamental physics and astrophysics are 

>0 ' being opened to scrutiny. 

I will also outline some of the ideas I advanced with Wheeler in our book "Black 
Holes Gravitational Waves and Cosmology"^ and in our article "Introducing the 
Black Hole'^(see Figs.[Tl[2l|3l|4]): they are coming to full fruition and new additional 
topics need urgent attention. 

It has been traditionally established that the understanding of an astrophysical 
phenomenon is only reached when its energetic aspects have been properly identi- 
fied. This has certainly occurred in three cases of clear success in the understanding 
of: 



(1) the stellar evolution, which was reached after thermonuclear reactions were fully 
tested in Earth-bound experiments and their results applied to astrophysics (see 
Refs.IMl; 

(2) pulsars, reached as soon as the rotational energy of neutron stars was recognized 
as the energy source of these phenomena (see Refs. fTSHTT]) : 

(3) the nature of binary X-ray sources, achieved as soon as the role of the gravita- 
tional energy in the fully general relativistic regime was established; the already 
well known results for a Schwarzschild metric were extended to the Kerr metric 
(see Fig. IH Ref. [IJ also quoted in Ref. El and Refs. I18|19p : these results were 
applied to the study of neutron stars and black holes in binary X-ray sources 
(seeRef.EQl). 
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Introducing the black hole 


Introducing the black hole 


^ 
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FIG. 2. (Reproduced from Ruffini and Wheeler, Ref. 4, with their kind permission.) Decay of a particle of rest- 
plus-ldnetic energy ^^ into a particle which is captured into the black hole with positive energy as judged locally, 
but negative energy J?i as judged from infinity, together with a particle of rest-plus-Idnetic energy E^ >E^ which 
escapes to infinity. The croas-hatched curves give the effective potential (gravitational plus centrifugal) defined 
by the solution B of Eq. (2) for constant values oip^ and i*. 

Fig. 1. Above: It was in the writing of our ESRO reporf-J with Johnny, which later became 
the first chapters of our book on "Black Holes, Gravitational Waves and Cosmology'ii' that some 
of the crucial new ideas on the last phases of gravitational collapse leading to a black hole were 
formulated. They were summarized in the celebrated article "Introducing the Black Hole" .121 There, 
the crucial new concept of a "black hole" , as a physical and astrophysical system and not just as an 
anal3rtic solution of the Einstein equations, was presented for the first time. Below: The particle 
decay process in the field of a black hole: figure and original caption reproduced from Ref. |4| 



If we turn from these very well understood and long duration phenomena to 
the most energetic and transient ones such as SNe and GRBs, it is well accepted 
that gravitational energy in the general relativistic regime plays the principal role 
in their explanation, namely the energy originating in the process of gravitational 
collapse either to a neutron star or to a black hole. Such a release of gravitational 
energy in supernovae is necessary to trigger the release of thermonuclear energy 
processes as well as the copious emission of neutrinos leading to the formation of 
a neutron star. Similarly, the more extreme general relativistic conditions created 
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Er^asphere of a rotating bUck hole. The reBJon b«tw<:en the ^urf^ce gf infinite 
redshiH (outer) and the event horizon (inner), itcre ^hown in 9 cutaway view, is 
called trie "ergospherc." when a parTlcJe disintegrates in this region end one at 
ttte fragments falls into Ihe blacii hole, the other Iragment C3n escape to infjnity 
wltll rr^ore rest plus kinetic energy than the original particle. Figure 5 



Fig. 2. The concept of the "ergosphere" was introduced^ as well as the concept of extraction of 
rotational and electromagnetic energy from a black hole with their corresponding thcrmodynamical 
analogies, further formalized in Refs. HI51 

in the collapse to a black hole trigger the vacuum polarization quantum process 
with a vast production of electron-positron pairs leading to the creation of a GRB. 
Nevertheless, some open issues remain in both cases. 



1.1. Some open issues on SNe 

For the case of SNe there is no better prototype then the Crab Nebula see Fig. \5\ 

The preliminary understanding of the role of a neutron star in the process of 
gravitational collapse leading to a supernova and the fundamental role of supernovae 
in the creation of cosmic rays were presented in a short and epochal paper by Baade 
and Zwicky.'251 Following the work of George Gamow,'23 Oppenheimer and VolkofipS 
gave the first general relativistic computations for neutron stars and introduced the 
existence of a critical mass against gravitational collapse. The continued gravita- 
tional collapse, solely described by general relativity for a star with mass larger than 
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Fig. 3. Wheeler described another crucial property of black holes in a very provocative pictorial 
representation illustrating the concept of the black hole uniqueness, namely the unprecedented 
situation in physics that the final output of the process of gravitational collapse was only dependent 
on three parameters, mass, charge, and angular momentum, independent of the original details 
of the precursor astrophysical source. We shall see shortly how this uniqueness concept plays a 
fundamental role in GRBs. 



the critical mass, was then developed by Oppenheimer and Snyder.^ Ganiow and 
Schocnberg221 identified the fundamental role of the neutrino and antineutrino emis- 
sion in order to dissipate the enormous thermal energy developed in the early phases 
of gravitational collapse. The neutrino-antineutrino emission, through the URCA 
process, introduced the essential cooling needed for the formation of a neutron star. 
In a pioneering work, Hoyle and Fowler^ described the more complex thermonu- 
clear process occurring in the gravitational collapse to a neutron star. They clearly 
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Fig. 4. Above: The effective potential experienced by a test particle moving in the equatorial 
plane of an extreme Kerr black hole. For corotating orbits, with positive values of the angular 
momentum, the maximum binding of 42.35% of the rest mass of the test particle is reached at the 
horizon. Details in Refs. fLlGl Below: Nuclear vs. gravitational binding energy in a Schwarzschild 
black hole compared and contrasted. The gravitational binding energy in the Kerr case is even 
bigger. See Ref. [T] also quoted in Ref. (6] 



differentiate two different possibilities: the gravitational collapse starting from a 
white dwarf, with a mass larger than the Stoner-Chandrasekhar critical mass}22l 
leading to a type I SN; or, alternatively, the gravitational collapse starting from the 
larger core of a more massive star, leading to a type II SN. The very extensive work 
on thermonuclear processes preceding and occurring at the onset of gravitational 
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Fig. 5. The Crab Nebula has certainly been one of the most complex systems ever analyzed by the 
human mind. Introduced by Messier in his famous catalog of nebulae^ as Ml, it was discovered 
by John Bevis and the name "Crab Nebula" was attributed by William Parsons, 3rd Earl of 
Rosse. The first images of its filamentary structure were obtained by Baade.l^ It was thanks to 
the fortunate interaction between Jaan Oorl|2Sl and the sinologist J.J. DuyvendaU23] that it was 
possible to identify the Crab nebula as the remnant of the SN which exploded in 1054 A.D. This 
event was recorded in the Chinese chronicles of the Sung dynasty and in the Japanese records 
Mei-Getsuki. This fascinating history has been reconstructed in the book by I.S. Shklovsky on 
SNeBU 



collapse has been exhaustively summarized in the classic book of David Arnett.'SS 
The observations of the supernova 1987 A together with the associated neutrinos 
was a particularly important step in the verification of this very complex scenario.'22 
There are nevertheless still open questions about the release of the gravitational 
energy: still unexplained is the process of the expulsion of the SN remnant. We still 
lack an understanding of such a fundamental process. The possibility of explaining 
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such an expulsion process by the very large flux of neutrinos emitted during the 
cooling of the collapsing core has been advocated by Jim Wilson and his coUabora- 
torsPSHSZl But it is a matter of fact that still today the theory is unable to explain 
this expulsion process.'23 

My collaborators and I, working in the field of relativistic astrophysics, are 
currently focused on seeing whether the way out of this impasse might be found 
in some crucial properties of neutron stars which have been neglected so far in a 
oversimplified approach. This crucial idea may well also have been missed in the 
extremely complex and time consuming numerical simulations. 

There is no doubt that there is plenty of gravitational energy available in the 
process of gravitational collapse. The gravitational energy can be estimated simply 
by (see e.g. Ref . l38l and references therein): 

where M and R are the mass and radius of the neutron star. If one assumes for the 
final radius of the collapsing star 

GM 

then the gravitational energy is 

E - -Mc^ . (3) 

a 

For 10 < a < 10^, typical for a neutron star, one readily obtains gravitational 
energies on the order of 10^^-10^'^ erg, which can well explain the observed energy 
of SNe in the range 10'*^-10^^ erg. In order to utilize this gravitational energy, either 
for the expulsion of the remnant or for the observed electromagnetic or particle 
energy flux emitted during the process of gravitational collapse, it appears essential 
to minimize the kinetic energy of the implosion of the core and transform the 
gravitational energy into an alternative form of energy with the necessary explosive 
power. There are two obvious mechanisms for obtaining this minimization process 
for the kinetic energy of the implosion; either the presence of rotation or the presence 
of new electrodynamical processes. The role of rotational energy is clearly revealed 
by the rotational energy of the newly formed neutron star, which is then released 
on time scales of thousands to millions of years in pulsars. 

Turning to alternative instantaneous braking mechanisms occurring during the 
process of gravitational collapse, we are currently exploring the possibility of using 
electrodynamical processes. In our analysis, even the electrodynamics of the neutron 
star configurations of equilibrium, especially in the interface of the crust to the core 
of the neutron star, needs further scrutiny. Until now the condition of local charge 
neutrality in the neutron star configuration of equilibrium has been adopted mainly 
for mathematical simplicity (see Sec. 14. 2|) . The current stringent energy release 
requirements, both in SN events and in GRBs, as we will see shortly, demand a 
deeper analysis of the electrodynamical properties during the gravitational collapse 
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Fig. 6. Left: the Vela 5A and 5B satellites. Right: a typical event as recorded by three of the 
Vela satellites. Details in Ref. 139 1 by I. Strong. 



process. A good starting point is the study of the equihbrium of a neutron star 
imposing global but not local charge neutrality. In such cores overcritical electric 
fields can develop and a large amount of electron-positron plasma can in principle 
be created during the process of their gravitational collapse, leading to explosive 
effects. This field of research will be summarized bricfiy in the conclusions. 



1.2. GRBs: The first steps of the fireshell vs. fireball scenario and 
the three paradigms 

For introducing the topic of GRBs there is no better illustration of this momentous 
research then the images of the original Vela Satellites (see Fig. |6|). The discovery 
of GRBs was a fortunate outcome of one of the most alarming episodes of the cold 
war between the United States of America and the then Soviet Union. The mis- 
sion was inspired by an unconventional proposal by our colleague Yacob Borisovich 
Zel'dovich.SSHlal xhe first public announcement of the GRB discovery occurred in 
a session I organized with Herb Gursky at the AAAS meeting in San Francisco.!^ 
On the meager observational evidence of the results presented in San Francisco, I 
worked out a model for GRBs with Thibault Damour, based on the possibility of 
explaining the origin of their gamma and X-ray fluxes by the process of vacuum 
polarization "a la Heisenberg-Euler-Schwinger" occurring around a black hole en- 
dowed with electromagnetic structure.'^ In that paper, written a few months after 
the announcement of the GRB discovery, we clearly pointed out three aspects: 
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(1) The role of the extractable electromagnetic energy of the black hole, or, more 
simply, the "blackholic" energy,Eil as the GRB energy source, implicit in the 
mass-energy formula of the black hole given by Christodoulou and RufflniP 






S = 167rm^,. , 



(4) 



^ + ifit,4 < 1 



4m^^ 16mf, 

^ SS = S2Trmir5mir > , 



where m is the total mass-energy of the black hole, rriir the irreducible mass, S 
the surface area and e and L are the black hole charge and angular momentum 
in geometrical units {G = c = 1). From this formula it follows that the amount 
of energy stored in a black hole, and in principle extractable by reversible trans- 
formations during the process of gravitational collapse, could be as high as 29% 
of the rotational energy and 50% of the electromagnetic energy of the black 
hole. The existence of what became known later as the "blackholic energy" 
extraction process (see Ref. l4Tj) had its beginnings here. In order to have an in- 
stantaneous explosive process, we directed our attention to the electromagnetic 
energy component in the mass-energy formula. We expected that the extraction 
of rotational energy, the other form of the blackholic energy, would take place 
on a much longer time scale, on the order of millions of years, e.g. in active 
galactic nuclei. 

(2) The creation of an electron-positron pair plasma by the vacuum polarization 
process in the field of a Kerr- Newman black hole, representing the actual process 
of energy extraction from the black hole; these processes of e+e" pair creation 
do approach the reversible transformations of a black hole introduced in Ref. [51 

(3) The identification of a typical "energy scale" on the order of ~ 10^'^Mbh/Mq 
erg, where Mbh is the black hole mass, as the typical energy of a GRB to be 
expected in our model. 

The theoretical background of this work was well grounded: 

• on the mathematical side the solution of the Einstein-Maxwell equations found 
by Roy Kerr,E51 and by Ted Newman and collaborators ,H3 as well as, on the 
physical side, the mass-energy formula of the Kerr-Newman black holes (see 
Eq.HD; ^ 

• on the quantum physics side, the Heisenberg-Euler'*^ and Schwingei^S formalism 
for the creation of an e~^e~ pair plasma by the vacuum polarization process; 
the basic equations in general relativity have been properly identified in Ref. |44] 
and fully reviewed in Ref. [49l 

• on the observational side, after the discovery by the Vela satellites a large mo- 
mentum was gained in the international community making GRBs a main object 
of the observations of a large variety of dedicated missions worldwide. One of 
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Fig. 7. The Coinpton Gamma-Ray Observatory satellite and the position in the sky of the ob- 
served GRBs in galactic coordinate. Different colors correspond to different intensities at the 
detector. There is almost perfect isotropy, both in the spatial and in the energetic distributions. 



the major issues addressed at the time was to know the actual location and 
distance of these sources in the universe. Just knowing their flux on the Earth, 
in the absence of such a distance determination, it was not possible to know 
their energetics. The GRB energies, in fact, could vary over an enormous range 
of values depending on their location in the solar system, in the galaxy, or in 
the distant extragalactic universe. Consequently, an enormous number of mod- 
els existed.^ This problem was finally overcome by the epochal observations by 
the BeppoSAX satellite, unequivocally confirming the typical "energy scale" I 
had previously introduced with T. Damour.'^ 

The Compton Gamma-Ray Observatory satellite, with the eight BATSE detec- 
tors, was the first of many very successful missions devoted to the study of GRBs. 
The outcome of these early observations led to two remarkable discoveries: the sub- 
stantial isotropy of the distribution of the sources in the sky (see Fig. [7]), as well as 
a clear separation of two families of events: the short GRBs, with observed duration 
shorter than ^ 1 seconds, and the long GRBs, corresponding to the complementary 
case (see Fig. [H). 

The number of dedicated missions soon increased (see Fig. |9]). The crucial dis- 
covery came from the Italian-Dutch satellite BeppoSAX, ^^ jointly operating, for 
the first time, in both the X-ray and gamma-ray domain. This allowed the "prompt 
radiation" observed in gamma-rays to be associated with a long-lasting (10^-10^ 
s) "afterglow" observed in X-rays. This extremely important discovery has led to 
a better understanding of the nature of the GRB sources and offered the essential 
tool for improving their localization in the sky. It triggered a coordinated and com- 
plementary set of observations in the optical and radio energy bands: in turn, the 
optical identification, made possible by the timely development of the most power- 
ful family of optical telescopes at Mauna Kea and at the VLT in Chile, allowed the 
determination of the GRB distance and therefore of their energetics (see Fig. |9]). 
This led to distances of the sources with z in the range between z ~ 0.0084 and 



July 6, 2011 0:9 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



11 



100 



OCT 

X 



m 
m 
<D 
C 

-a 
u 

T3 
(U 
OB 

to 

> 
CO 
I 

0) 

o 
c 

cu 
3 



oo 
c 



10 



0.1 



n I I I I 1 1 



m]- 




_i I 



_i 



0.01 



0.1 



10 

(sec) 



100 



1000 



Fig. 8. The energy fluence- averaged hardness ratio for short (T < 1) and long (T > 1 s) GRBs 
are represented. Reproduced with the kind permission of M. Tavani, from Ref. [52] where the details 
are given. 



z ~ 8 2p2t^ and corresponding values of the energetics all the way up to the order 
of 10^'* erg, exactly in the range we had predicted with Thibault Damour 25 years 
earlier. 

We soon returned to our model by introducing three major new conceptual 
developments: 

(1) The identification of the region around the black hole where an overcritical elec- 
tric field could occur, leading to the creation of e^e^ pairs out of the vacuum. 



The 






were assumed to be in thermal equilibrium; the first identification 



was performed around a Reissner-Nordstr0m black hole. From the Greek name 
"5va<;, 5vaSo<;^^ for "pairs" , the name "dyadosphcre" was introducecP^I^ for 
this region, characterized by a total energy E'fot, a mean energy of the pair 
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Fig. 9. The large flotilla of space observatories which followed the Vela satellites is here repre- 
sented by the sequence of the Compton satellite and the fundamental contribution by the Bep- 
poSAX satellite which, through its discovery of the afterglow, has allowed the determination of a 
sharper position for GRBs and consequently the possibility to identify their optical and radio coun- 
terparts with many space and Earth-based telescopes. This has allowed the determination of their 
cosmological distance and of their isotropic equivalent energy on the order of lO^** erg. Currently, 
the satellites XMM, Chandra, Swift, Fermi and Agile are giving new fundamental contributions 
on the energy, time variability and spectra of GRBs. 



created on the order of 1-4 MeV and a characteristic radius 10*-10^ cm. For 
details see Ref. [ST] 

(2) The study with Jim Wilson of the dynamics of the e~^e~ pairs as the fun- 
damental acceleration process of GRBs, reaching bulk Lorentz gamma factors 
7 ~ 10^-10^ (the pair-electromagnetic pulse, PEM pulse). This theoretical re- 
sult was obtained in Ref. EH] (for details see Sec. [3]). 

(3) Still with Jim Wilson, the study of the acceleration process of the optically thick 
e^e^-baryon plasma (the pair-electromagnetic-baryon pulse, PEMB pulse); this 
theoretical result was published in Ref. [591 This analysis introduced the new 
concept of plasma baryon loading defined by i? = MBC^/E^^f, where Mb is the 
total mass of the baryons. To our surprise, we found a maximum value of the 
baryon loading for the existence of the electron-positron acceleration process 
B < 10^2 (for details see Sec.E]). 
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Some of these considerations have been consistent with and others in contrast 
to other articles in the contemporary hterature: 

(1) The role of the e^e" plasma, originally presented in our work with Damourp2 
was soon generally adopted, with and without references Ip^lM] 

(2) There has been a substantial difference between our approach, purporting an 
initial thermal distribution of the e+e~ plasma in the dyadosphere, and the 
work by Cavallo and Rees,!^ purporting the total annihilation of the e+e" 
pairs originating in a process of gravitational collapse; 

(3) The model "a la" Cavallo-Rees has led to the concept of a "fireball," a very 
hot cavity originating from the annihilation of the e^e~ pairs pushing on the 
surrounding circumstellar matter. In contrast, our model has led to the alter- 
native concept of the "fireshell" in which the gradual annihilation of the e+e" 
pairs lead to a self-acceleration of an optically thick shelliSSEH xhe cavity inside 
such a shell is practically at zero temperature, and the e+e~ pairs gradually 
annihilate lasting until the fireshell becomes optically thiniSS 

To make this general situation even more exciting, a totally unexpected phe- 
nomena has been observed: the discovery of a GRB family that is relatively weak 
(10''^-10^^ erg) and close {z < 0.17), coincident in space and time with a SN event 
(see e.g. Ref . l67l and references therein). 

The "majority" point of view reached consensus on three main points: 

(I) That short GRBs originate from binary mergers of white dwarfs, neutron stars 

and/or black holes in all possible combinations .1S2M1 
(II) That long GRBs originate from the collapse of massive stars and the observed 
spikes in their light curves originate in the prolonged activity of an inner 
engine .'S2HI2] jjj order to make the energy requirement less stringent they pos- 
tulated the existence of a jet structurep^Ell jn particular, the existence of a 
standard energetics for all GRBs of ~ 10^^ erg^ was also claimed. 
(Ill) The observations of SNe associated only with long GRBs was considered a 
clear support for the idea that all long GRBs originate from supernovae in the 
collapse of very massive stars: the "collapsar" modeL'^^'^ 

Some of aspects of our work have offered a different point of view. We expressed 
them in three basic paradigms: 

(1) The relative space-time transformation (RSTT) paradigmpS the first paradigm 
emphasizes the relevance of having the correct space-time parametrization of 
the source, which necessarily implies the knowledge of the equations of motion 
of the system and its entire worldline. This procedure was in contrast with the 
current practice of describing the GRB nature by a piecewise analysis. 

(2) The interpretation of the burst structure (IBS) paradigm!^ the second 
paradigm emphasizes the structure of the canonical GRB as composed of a 
proper GRB (P-GRB) and an extended afterglow, whose characteristics are 
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mainly dominated by the two parameters of the total energy Etot and of the 
baryon loading B — M bc^ / Etot of the e+e^ plasma (where AIb is the mass of 
the baryon loading): in the limit of S — > the canonical GRB would result in 
a short GRB and, in the opposite limit of B — >■ 10^^, which is the maximum 
possible valuepS it would result in a long GRB. Our model assumes spherical 
symmetry. 
(3) The GRB-supernova time sequence (GSTS) paradigmfSSl the third paradigm 
introduces the concept of "induced gravitational collapse," which considered 
the effect of a GRB triggering the late phase evolution of a highly evolved 
companion star leading to a supernova (SN). 

As we will see, the consistent application of these paradigms has led to con- 
clusions quite different from the ones (I), (II), (III) formulated by the majority, 
especially with reference to short GRBs and to the collapsar model. 

Before introducing the recent developments in our model, I outline some general 
considerations out of first principles establishing an upper limit on the black hole 
mass and, consequently, on the energetics of GRBs. I also recall a scenario outlined 
with John Wheeler. 

2. General considerations on GRB energetics 

The observations of GRBs, made possible by the large flotilla (see Fig. [9]) of space 
observatories in the gamma-ray, X-ray and visible energy ranges as well as by Earth- 
bound telescopes in the visible and radio energies, has led to the discovery of a GRB 
almost every day. We exemplify in Fig. [10] some recent representative observations 
of GRBs, with their energetics falling into the range 10'*^-10^^ erg. 

2.1. The upper limit to GRB energetics 

It is appropriate to return to an order of magnitude estimate to express the ne- 
cessity of having a black hole of less than ^ IQMq in order to explain the ener- 
getics of GRBs. From the expression for the gravitational energy in the process of 
gravitational collapse, given by Eq. ([T]), and assuming that the dimension of the 
gravitationally collapsed object is on the order of i? = aGM/c^ , where a is on the 
order of 1 or 2, typical of a black hole, the gravitational energy available is on the 
order of Mf? j a. This simple order of magnitude estimate is in agreement with the 
rigorous energetics defined by the mass-energy formula of the black hole (see Eq. S]). 
If one takes the average density of the collapsing object as: 

Pec = A//[(4/3)^a3(c'/c2)3M3] , (5) 

and expresses it in terms of the mass M* = Afi/m^, one obtains: 

Pec = TO„/[(4/3)7ra3(M/M*)2(n3/c3m3)] . (6) 

If one takes the nuclear density 

p„ = Am„/[(4/3)^(nVc^™')] , (7) 
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Fig. 10. A selected number of GRBs with energy in the range 10*^-10^^ erg are represented on 
our past hght cone as a function of their cosinological redshift and of the time from the big bang. 
The cosmological horizon is represented by the dashed curve. 

in order to have a vacuum polarization process occur, we must have pcc > Pn- This 
gives an upper hmit to a black hole mass of: 



Mbh < {mn/m^f^^/{M*VA^) - l/Aa^/2 - 9M( 



7 



(8) 



The collapsing core should have a mass on the order of IOMq. Correspondingly 
there exists an absolute upper limit to the GRB energetics at Eqrb ^ 10^^ erg. 



2.2. The role of the Kerr-Newman solution and open issues in 
gravitational collapse 

The great advantage of the existence of the Kerr-Newman metric is that we have 
an exact mathematical solution of the Einstein-Maxwell equations in which we can 
probe the order of magnitude estimate of the amount of energy released, the spectral 
distribution of the radiation, and the observed characteristic time of the process at 
infinity to be expected in a gravitational collapse. It is clear, however, that the use of 
a Kerr-Newman solution is made just for mathematical and idealized physical con- 
venience. The real GRB astrophysical description will not have an already formed 
black hole. It will need, in reality, the detailed description of the process of gravita- 
tional collapse to a black hole: possibly the most difficult process to be described in 
physics and astrophysics. For understanding the black holes as energy sources it is 
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essential to identify the transient process leading to building up the electromagnetic 
structure of the black hole forming process during the gravitational collapse. The 
build-up process of an overcritical electric field is ended in a sudden discharge by 
the quantum process leading to the formation of e~^e~ pairs I^^M] 'When it comes 
to the detailed description of such a phenomenon, the necessity of broadening our 
current knowledge to new theoretical fields appears to be necessary. Even the sim- 
plest concept of a neutron star has to be generalized by properly accounting for new 
electrodynamical processes^ (see Sec. l4.2p . To analyze their dynamical evolution in 
the process of gravitational collapse is still a formidable theoretical challenge due 
to the interconnections of general relativity and all the relativistic field theories. 
The exploration of this much more complex physical and astrophysical reality and 
the understanding of this "terra incognita" is certainly one of the most exciting 
fields of current research. What we can certainly ascertain since now is that we are 
obtaining a more unified physical picture extending from the microphysical world of 
nuclear physics all the way to general relativistic astrophysical systems taking into 
due account all fundamental interactions. This novel scientific experience will likely 
influence also a new approach to cosmology as soon as the cosmological density 
overcomes nuclear densities. 

On the other hand the formation of the horizon of the final almost Schwarzschild 
black hole, occurring after the annihilation of the electromagnetic structure of the 
Kerr-Newman solution into the e+e^ plasma creation, is certainly mathematically 
interesting but physically of smaller relevance since, by its own definition, no black- 
holic energy is available and it is therefore only indirectly observable by some ac- 
cretion phenomena. 

Progress toward the goal of understanding the dynamics of gravitational col- 
lapse is being made (see e.g. Ref . H^ . The basic relevance of the Coulomb structure 
in reducing the kinetic energy of gravitational collapse and in creating an electro- 
dynamical structure around the black hole has been highlighted in a set of papers 
using simplified models (see Fig. [TT]as well as Ref. W2\ and references therein). 

The general scenario can be summarized as follows: we start with an initial 
configuration having global but not local charge neutrality. In the process of gravi- 
tational collapse of such a configuration there is an increase of the electrodynamical 
structure. The formation of overcritical electric fields leads to the electron-positron 
plasma creation and the disappearance of all the electromagnetic structure.*^ Fi- 
nally, a neutral black hole is left over by the process of gravitational collapse. An 
estimate on the order of magnitude of this process has been simulated by con- 
sidering the decay of an already formed Kerr-Newman black hole endowed with 
an overcritical electric field, 1^ to an almost Schwarzschild black hole without any 
electromagnetic structure. 
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Fig. 11. Energetics of a shell such that Afi = Mq, for selected values of the charge. In the first 
diagram Q = 0; the dotted-dashed line represents the total energy for a gravitational collapse 
without any braking process as a function of the radius R of the shell; the solid step function line 
represents a collapse with suitable braking of the kinetic energy of implosion at selected radii; the 
dotted line represents the rest mass energy plus the gravitational binding energy. In the second and 
third diagram Q/Mq = 0.7, Q/Mq = 1 respectively; the dotted-dashed and the dotted lines have 
the same meaning as above; the solid lines represent the total energy minus the kinetic energy. 
The region between the solid line and the dotted line corresponds to the stored electromagnetic 
energy. The region between the dotted-dashed line and the solid line corresponds to the kinetic 
energy of collapse. In all the cases the sum of the kinetic energy and the electromagnetic energy 
at the horizon is 50% of Mq. Both the electromagnetic and the kinetic energy are in principle 
extractable. It is most remarkable that the same underlying process occurs in the three cases: the 
role of the electromagnetic interaction is twofold: a) to reduce the kinetic energy of implosion by 
the Coulomb repulsion of the shell; b) to store such an energy in the region around the black hole. 
The stored electromagnetic energy is extractable as shown in Ref. l82l and leads to the pair creation 
process if the field is overcritical. See also Ref. 1831 



2.3. The black hole uniqueness theorem in GRBs 

In Fig. [3] we have recalled the exotic way of illustrating the black hole uniqueness 
theorem advanced by John Wheeler and expressed in his characteristic pictorial Ian- 
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Fig. 12. The black hole uniqueness theorem represented by John Wheeler in Fig.[3]is here applied 
to the case of GRB emission. 



guage. Indeed, GRBs offer a most grandiose scenario to check and test the validity of 
this uniqueness theorem. As we have mentioned above, for reasons of energetics, the 
GRB energy emission process occurs in the latest phases of gravitational collapse, 
when the black hole horizon is approached. The GRB progenitors can originate in a 
vast range of possibilities: from X-ray binaries, binary neutron stars, single star col- 
lapse, black holes in the core of globular clusters, or intermediate mass black holes. 
As recalled above, the process of emission by GRBs occurs in the latest phases of 
gravitational collapse, when the black hole horizon is approached: there the major- 
ity of energy emission occurs, and this emission will be quite independent of the 
nature of the source. It will be uniquely characterized by the mass, the angular 
momentum and the electrodynamical structure of the collapsing core. In this sense, 
we can afford to make a unique model for all GRBs. The observed GRB properties 
will be only different in the time scale and spectral and energetic properties of the 
black hole "dyadotorus" , which will be defined shortly (see Fig. [T2|) . The effective 
electrodynamical structure of the collapsing core, idealized in our model as an al- 
ready formed black hole, will generally depend on a self-consistent treatment of the 
strong, weak, electromagnetic and gravitational interaction of the collapsing core, 
and will also depend on the fermion statistical properties (see Sec. l4.2|) . The essen- 
tial information on this collapsing phase will be encoded in the structure, spectrum 
and time evolution of the P-GRB or of the short GRBs, see Sec. [H 
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Black Holes, Gravitational Waves and Cosmology 




Figure 33 A rotating star with dense core A collapses to a pancake neutron 
star B; it fragments C; the fragments lose energy in periodic and splash 
gravitational radiation and recombine. The lower curve gives a schematic 
representation of the quadrupole moment as a function of time. Between B 
and C impulse radiation is created in the act of fragmentation not adequately 
described by the one indicated component of the quadrupole moment 
tensor. Between C and D multiply periodic radiation is given out until 
at D two fragments have lost enough angular momentum so that they 
combine with a splash of gravitational radiation ; similarly at E, etc. 

Fig. 13. Reproduced from Ref. 'S' See also on the self-gravitating rotating configurations Ref. 1841 
in collaboration with Simonetta Filippi and references therein. 



I finally recall the pursuit and plunge scenario (see Fig. I13p. which may become 
relevant in the study of possible multiple GRB phenomena. 
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2.4. Summary of next sections 

I am going to outline in the next sections some of the most significant moments in 
our research: 



(1) the "canonical" GRB predicted by the fireshell model; 

(2) some basic progress in the theoretical understanding of fundamental physical 
processes motivated by the study of GRBs; 

(3) a new class of GRBs identified thanks to the application of the above-mentioned 
three paradigms. 



3. The canonical GRB scenario in the fireshell model 

We assume, within the fireshell model, that all GRBs originate from an optically 
thick e+e" plasma with total energy Ef^^. in the range lO^'^-lO^"' erg, a temperature 
T in the range 1-4 MeV and typical radii 10^-10^ cml^Such an e+e" plasma has 
been widely adopted in the current literature (see e.g. Refs. I68I87I and references 
therein) . 

The dynamics of the GRB develops in two very different phases: the first one, 
optically thick, and the second one, optically thin. 

In the initial optically thick phase, after an early expansion, the e^e^-photon 
plasma reaches thermal equilibrium with the engulfed baryonic matter Mb de- 
scribed by the dimensionless parameter B — Mbc^ /Ef^f. that must be in the 
range B < iQ-^pSMi xhe fireshell composed of e+e^-photon-baryon plasma self- 
accelerates to ultrarelativistic velocities and finally reaches the transparency con- 
dition. A flash of radiation is then emitted. This is the P- GRB .^^Sl Different current 
theoretical treatments of these early expansion phases of GRBs are compared and 
contrasted in Refs . [88l and [89l The amount of energy radiated in the P-GRB is only 
a fraction of the initial energy Ef^^ . 

At the decoupling between matter and radiation at the transparency point, 
the optically thin phase starts, characterized by a baryonic and leptonic matter 
fireshell endowed with ultrarelativistic Lorentz gamma factors (200 < 7 < 2000). 
Such an accelerated optically thin fireshell gives rise to a multi-wavelength emission 
by inelastic collisions with the circumburst medium (CBM). This is the extended 
afterglow. It has three different regimes: a rising part, a peak and a decaying tail. 
We therefore define a "canonical GRB" light curve with two sharply different com- 
ponents (see Fig.[14l):'^9 ^^ ^^'^'^ 1) the P-GRB and 2) the extended afterglow. What 
is usually called "prompt emission" in the current literature, in our canonical GRB 
scenario is therefore composed of the P-GRB together with the rising part and the 
peak of the extended afterglow. The unjustified mixing of these two components, 
which originates from different physical processes, leads to difficulties in the current 
models of GRBs. 
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Fig. 14. The "canonical GRB" light curve computed theoretically for GRB 991216. The prompt 
emission observed by BATSE is identified with the peak of the extended afterglow, while the 
small precursor is identified with the P-GRB. For this source we have _E'J* = 4.83 X 10^^ erg, 
S ~ 3.0 X 10~^ and {n^t,n) ~ 1-0 particles/cm^. Details in Refs. I79l85l86l 



3.1. The optically thick phase 

In Fig. [TS] we recall the evolution of the optically thick fireshell Lorentz gamma 
factor as a function of the external radius for 7 different values of the fireshell 
baryon loading B and two selected limiting values of the total energy ii^*°* of the 
e^e^ plasma. We can clearly identify three different eras: 
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(1) Era I: The fircshell is made up only of electrons, positrons and photons in 
thermodynamic equilibrium (the "pair-electromagnetic pulse" , or PEM pulse 
for short). It self-accelerates and begins its expansion into vacuum, because the 
environment has been cleared by the black hole collapse. The Lorentz gamma 
factor increases with radius and the dynamics can be described by the energy 
conservation and the condition of adiabatic expansion.ESMl 

(2) Era II: The fircshell impacts with the non-collapsed bayonic remnants and en- 
gulfs them. The Lorentz gamma factor drops. The dynamics of this era can be 
described by imposing energy and momentum conservation during the fully in- 
elastic collision between the fircshell and the baryonic remnant. For the fireshell 
solution to be still valid, B < 10~^ must hold.'^ 

(3) Era III: The fireshell is now made up of electrons, positrons, baryons and pho- 
tons in thermodynamic equilibrium (the "pair-electromagnetic-baryonic pulse" , 
or PEMB pulse for short). It self-accelerates again and the Lorentz gamma 
factor increases again with radius up until when the transparency condition 
is reached, going to an asymptotic value jasym = l/B. If B ^ 10^^ the 
transparency condition is reached when 7 ~ "fasym- On the other hand, when 
B < 10~^, the transparency condition is reached much before 7 reaches its 
asymptotic value. ^^"^^ In this era the contribution of the rate equation starts to 
be crucial in describing the annihilation of the e+e^ pairs: 

B ^ f)v ^ 

|^e± = -^e± y ^ + av-, {n!^ (T) - N!^) , (9) 

where 7Vg± is the number of e~^e~ pairs and 7Vg±(r) is the number of e^e~ 
pairs at thermal equilibrium at temperature yElMl 

In the "fireball" model in the current literature the baryons are usually considered 
to be present in the plasma from the very beginning. In other words, in the fire- 
ball dynamics there is only one era corresponding to the Era III abovc^^ilMIMIMl 
Moreover, the rate equation is usually neglected, and this affects the reaching of the 
transparency condition. A detailed comparison between the different approaches is 
given in Ref . [88l 



3.2. The transparency point 

At the transparency point, the value of the B parameter governs the ratio between 
the energetics of the P-GRB and the kinetic energy of the baryonic and leptonic com- 
ponents giving rise to the extended afterglow. It governs as well the time separation 
between the peak luminosities of the P-GRB and of the extended afterglow.'^^IMl 

By solving the rate equation we have evaluated the evolution of the temperature 
during the fireshell expansion, all the way up to when the transparency condition 
is reached.'^^m 

In the upper left panel of Fig. [16] we plot, as a function of B, the fireshell radius 
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Fig. 15. The Lorentz gamma factor of the expanding firesheli is plotted as a function of its 
external radius for 7 different values of the firesheli baryon loading B, ranging from B = 10~* 
and B = 10"^, and two selected limiting values of the total energy i?*j* of the e+e~ plasma: 
E*-'^ = 1.17 X 10*^ erg (upper panel) and E^'^ = 1.47 X 10^'' erg (lower panel). The asymptotic 
values 7 — > 1/B are also plotted (dashed horizontal lines). The lines are plotted up to when the 
firesheli transparency is reached. For details see Ref. 1591 



at the transparency point. The plot is drawn for four different values of £'*±* in the 
interval [10'*^, 10^^] erg which encompasses the observed isotropic GRB energies. 
In the upper right panel of Fig. [T6]we plot, as a function of B, the firesheli tern- 
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Fig. 16. At the fireshell transparency point, for 4 different values of i?*°'i ^^ plot ^^ ^ function of 
B: (upper left) the fireshell radius in the laboratory frame; (upper right) the fireshell temperature 
in the co-moving frame T^°'" (thicker lines) and the one Doppler blue-shifted along the line of 
sight toward the observer in the source cosmological rest frame T°^" (thinner lines); (lower left) 
the fireshell Lorentz gamma factor 70 together with the asymptotic value 70 = l/S; (lower right) 
the energy radiated in the P-GRB (thinner lines, rising when B decreases) and the one converted 
into baryonic kinetic energy and later emitted in the extended afterglow (thicker lines, rising when 
B increases), in units of i?*"'. For details see Refs. [591791 



perature To at the transparency point, i.e. the temperature of the P-GRB radiation. 
The plot is drawn for the same four different values of i5*°* of the upper panel. We 
plot both the value in the co- moving frame J"™™ and the value Doppler blue-shifted 
toward the observed T"''* = (1 -I- /3o)7oT'o°™, where (3o is the fireshell speed at the 
transparency point in units of c\^ 

In the lower left panel of Fig.[Tn]we plot, as a function of B, the fireshell Lorentz 
gamma factor at the transparency point 70. The plot is drawn for the same four 
different values of £■*"* as in the upper panel. Also plotted is the asymptotic value 
7o = l/S, which corresponds to the condition when the entire initial internal energy 
of the plasma i?*°* has been converted into kinetic energy of the baryons.'SS! We see 
that such an asymptotic value is approached for B — > 10^^. We see also that, if 
E*J± increases, the maximum values of 70 are higher and they are reached for lower 
values of B. 

In the lower right panel of Fig. [16] we plot, as a function of B, the total energy 
radiated at the transparency point in the P-GRB and that converted into baryonic 
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Fig. 17. Here the energies emitted in the P-GRB (sohd hne) and in the extended afterglow (dashed 
Hne), in units of the total energy of the plasma, are plotted as functions of the B parameter for a 
typical value of _E*°' ~ 10^^ erg (see lower panel of Fig. llGt . When B < 10~^, the P-GRB becomes 
predominant over the extended afterglow, giving rise to a "genuine" short GRB. The figure also 
shows the values of the B parameter corresponding to some GRBs we analyzed, all belonging to 
the class of long GRBs. 



and leptonic kinetic energy and later emitted in the extended afterglow. The plot 
is drawn for the same four different values of i?*±* as in the upper panel and middle 
panels. We see that for B < 10~^ the total energy emitted in the P-GRB is always 
larger than that emitted in the extended afterglow. The limit i? — > gives rise to a 
"genuine" short GRB (see also Fig. [T7| , namely a GRB whose prompt emission is 
dominated by the P-GRB and which is followed by a very tiny extended afterglow, 
if any at all. On the other hand, for 3.0 x 10~^ ^ B < 10~^ the total energy emitted 
in the P-GRB is always smaller than the one emitted in the extended afterglow. 
If it is not below the instrumental threshold and if ricbm ~ 1 particle/cm'^, the 
P-GRB can be observed in this case as a small pulse preceding the main GRB 
event (which coincides with the peak of the extended afterglow), i.e. as a GRB 

"precursor" I79l89l90l92lt94l 

Particularly relevant for the new era of the Agile and Fermi satellites are the 
GRBs with B < 10"'^. In this case the P-GRB emission has an observed temperature 
up to 10'^ keV or higher. This high-energy emission has been unobservable by the 
Swift satellite. 

We must emphasize that all the above estimates have been done for a Reissner- 
Nordstrom black hole endowed with an overcritical electric field. Some differences 
should exist in the structure of the P-GRB if, instead of a dyadosphere, the col- 
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lapse gives birth to a dyadotorus. In addition, in the actual process of gravitational 
collapse the dyadosphere-dyadotorus formation will not occur in an asymptotically 
flat space. Its boundaries will be characterized by the characteristic processes of 
gravitational collapse. There is therefore the very exciting possibility that the ac- 
tual details of the process of gravitational collapse can be inferred in principle from 
the structure of the P-GRB or of the "genuine" short GRBs: from their duration, 
spectra and time variability. 

3.3. The optically thin phase 

The dynamics of the optically thin fireshell of baryonic matter propagating in the 
CBM has been obtained from the relativistic conservation laws of energy and mo- 
mentum (see e.g. Ref. 95). Such conservation laws are used both in our approach 
and in others in the current literature (see e.g. Refs. r72l86l95f[57|) . The main differ- 
ence is that in the current literature an ultra-relativistic approximation, following 
the Blandford-McKee^Sl self-similar solution, is widely adopted, leading to a sim- 
ple constant-index power-law relation between the Lorentz gamma factor of the 
optically thin "fireshell" and its radius: 

7 ex r-'^ , (10) 

with a = 3 in the fully radiative case and a — 3/2 in the adiabatic case.'^^EZl 
On the contrary, we use the exact solutions of the equations of motion of the 
fireshell.'^^ISMZMIiool ^ detailed comparison between the equations used in the 
two approaches has been presented in Refs. I95|97I99I100I In particular, Ref. |97] 
shows that the regime represented in Eq. ()10p is reached only asymptotically when 
7o ^ 7 ^ 1 in the fully radiative regime and 7^ » 7^ :^ 1 in the adiabatic regime, 
where 70 is the initial Lorentz gamma factor of the optically thin fireshell. 

In Fig. [TH] we show the differences between the two approaches. In the upper 
panel there are plotted the exact solutions for the fireshell dynamics in the fully ra- 
diative and adiabatic cases. In the lower panel we plot the corresponding "effective" 
power-law index a^ff , defined as the index of the power- law tangent to the exact 
solutionis 

fle// = --71 ■ (11) 

dlnr 
Such an "effective" power-law index of the exact solution smoothly varies from to 
a maximum value which is always smaller than 3 or 3/2, in the fully radiative and 
adiabatic cases respectively, and finally decreases back to (see Fig. [ 



3.4. Extended afterglow luminosity and spectra 

The extended afterglow luminosity in the different energy bands is governed by 
two quantities associated with the environment. Within the fireshell model, these 
are the effective CBM density profile, Ucbm, and the ratio between the effective 
emitting area Aejf and the total area Atot of the expanding baryonic and leptonic 
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Radial coordinate (r) (cm) 

Fig. 18. In the upper panel, the analytic behavior of the Lorentz gamma factor during the ex- 
tended afterglow era is plotted versus the radial coordinate of the expanding optically thin fireshell 
in the fully radiative case (solid line) and in the adiabatic case (dotted line) starting from 70 = 10^ 
and the same initial conditions as the GRB 991216.^23 In the lower panel are plotted the corre- 
sponding values of the "effective" power-law index Qe// (see Eq. mil l), which is clearly not constant 
but highly varying and systematically lower than the constant values 3 and 3/2 purported in the 
current literature (horizontal thin dotted lines). 



shell, 7?. = Aeff/Atot- This last parameter takes into account the CBM filamentary 
structur d^°^ l ^°^ l and the possible occurrence of fragmentation in the shell. 1^221 

Within the "fireshell" model, in addition to the determination of the baryon 
loading, it is therefore possible to infer a detailed description of the CBM, its aver- 
age density and its porosity and filamentary structure, all the way from the black 
hole horizon out to a distance r < 10^^ cm. Typical dimensions of these clouds of 
overdense material in the CBM are on the order of 10^^-10^^ cm and they show 
a density contrast on the order of Ancbm/incbm) ~ 10. A fascinating possibility 
for the origin of these clouds has been suggested by David Arnett ,'^22 namely that 
these clouds correspond to matter ejected in the latest phases of the thermonuclear 
evolution of the progenitor star, leading to the black hole formation originating 
the GRB. The interaction between the baryons and leptons of the relativistically 
expanding shell with these CBM clouds corresponds to the spikes in the gamma 
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between the P-GRB and the peak of the extended afterglow (i.e. the "quiescent time between the 
"precursor" and the main GRB event), measured in the source cosmological rest frame. The com- 
putation has been performed assuming a constant value of the CBM density in four different cases: 
^cbm = 1-0 particles/cm'^, n^,;,^ = 1.0 X 10~^ particles/cm^, n^j^ = 1.0 X 10~^ particles/cm^, 
"c6m = f-0 X 10~^ particles/cm^. The points represent the actual numerically computed values, 
connected by straight line segments. See details in Ref. 1421 



and X-ray light curve of the prompt emission. This description is missing in the 
traditional model based on synchrotron emission. In fact, the attempt to use the 
internal shock model for the prompt emission (see e.g. Refs. I65I68I87]) only applies 
to regions where r > 10^^ cm.^^SS] 

In our hypothesis, the emission from the baryonic and leptonic matter shell is 
spherically symmetric. This allows us to assume, in a first approximation, a modeling 
of the CBM distribution by thin spherical shells and consequently to consider just 
its radial dependence. l^S 

For simplicity and in order to have an estimate of the energetics, the emission 
process is postulated to be thermal in the co-moving frame of the shell.U^ The 
observed GRB nonthermal spectral shape is due to the convolution of an infinite 
number of thermal spectra with different temperatures and different Lorcntz and 
Doppler factors. Such a convolution is to be performed over the surfaces of constant 
arrival time of the photons at the detector (equitemporal surfaces, EQTSs; see e.g. 
Ref. [95|) encompassing the entire observation time interval .'^^Sl "We are currently 
considering additional effects as an alternative to a purely thermal spectrum.'^SIEnil 
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Fig. 20. We plot three theoretical extended afterglow bolometric light curves together with the 
corresponding P-GRB peak luminosities (the horizontal segments). The computations have been 
performed assuming the same B*^* and CBM structure as the GRB 991216 and three different 
values of B. The P-GRBs have been assumed to have the same duration in the three cases, i.e. 
5 s. For B decreasing, the extended afterglow light curve squeezes itself onto the P-GRB curve 
(Simulations by L. Caito, see details in Ref. I42|l . 



In Fig.[T9]we plot, as a function of S, the arrival time separation Ata between the 
P-GRB and the peak of the extended afterglow measured in the cosmological rest 
frame of the source. Such a time separation Ata is the "quiescent time" between the 
precursor (i.e. the P-GRB) and the main GRB event (i.e. the peak of the extended 
afterglow). The plot is drawn for the same four different values of i5*°* of Fig. [161 The 
arrival time of the peak of the extended afterglow emission depends on the detailed 
profile of the CBM density. In this plot a constant CBM density has been assumed 
in four different cases: ricbm = 1-0 particles/cm'^, ricbm = 1-0 x 10~^ particles/cm'^, 
ncbm = 1-0 X 10^*^ particles/cm^, Ucbm — 1-0 x 10~^ particles/cm^. We can see that, 
for 3.0 X 10^^ ^ B < 10^^, which is the condition for P-GRBs to be "precursors" 
(see above), Ata increases both with B and with i<^*°*- We can have Ata > 10^ s 
and, in some extreme cases even Ai^ '^ 10'^ s. For B < 3.0 x 10"'', instead, Ata 
shows a behavior which qualitatively follows the opposite of 70 (see the middle panel 
of Fig. [111). 

Finally, in Fig. [20l we present three theoretical extended afterglow bolometric 
light curves together with the corresponding P-GRB peak luminosities for three 
different values of B. The duration of the P-GRBs has been assumed to be the 
same in the three cases (i.e. 5 s). The computations have been performed assuming 
the same El^ and the same detailed CBM density profile as GRB 991216.'22! In 
this picture we clearly see how, for B decreasing, the extended afterglow light curve 
"squeezes" itself onto the P-GRB and the P-GRB peak luminosity increases. We 
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are currently trying to identify some GRBs having this feature. 

The radiation viewed in the co-moving frame of the accelerated baryonic matter 
is assumed for simplicity to have a thermal spectrum and to be produced by the in- 
teraction of the CBM with the front of the expanding baryonic shellP^ In Ref. 11061 
it was shown that, although the instantaneous spectrum in the co-moving frame of 
the optically thin fireshell is thermal, the shape of the final instantaneous spectrum 
in the laboratory frame is nonthermal. In fact, as explained in Ref. llOli the tem- 
perature of the fireshell is evolving with the co-moving time and, therefore, each 
single instantaneous spectrum is the result of an integration of hundreds of thermal 
spectra with different temperatures over the corresponding EQTS. This calculation 
produces a nonthermal instantaneous spectrum in the observer frame.'^^ 

Another distinguishing feature of the GRB spectra which is also explained 
within the fireshell model is their hard to soft transition during the evolution of 
the event jMloaKliU j^ f^^t the peak of the energy distributions Ep drift monoton- 
ically to softer frequencies with time.^^^ This feature explains the change in the 
power-law low energy spectral indejipJ^ a which at the beginning of the prompt 
emission of the burst (if = 2 s) is a = 0.75, and progressively decreases for later 
times.'^^ In this way the link between Ep and a identified in Ref. 11091 is explic- 
itly shown. This is due to the decrease of the Lorentz gamma factor and of the 
temperature in the co-moving frame .'^^EMl 

The time-integrated observed GRB spectra show a clear power-law behavior. 
Within a different framework (see e.g. Ref. 1113| and references therein) it has been 
argued that it is possible to obtain such a power-law spectra from a convolution 
of many non-power-law instantaneous spectra monotonically evolving in time. This 
result was recalled and applied to GRBsIl^ assuming for the instantaneous spectra a 
thermal shape with a time varying temperature. It was shown that the integration 
of such energy distributions over the observation time gives a typical power-law 
shape consistent with the GRB spectra. 

Our specific quantitative model is more complicated than the one considered in 
Ref. 11141 the instantaneous spectrum in the Fireshell model is not a black body. 
Each instantaneous spectrum is obtained by an integration over the corresponding 
EQTS ,1251221 themselves a convolution, weighted by appropriate Lorentz and Doppler 
factors, of ~ 10^ thermal spectra with time varying temperature. Therefore, the 
time-integrated spectra are not plain convolutions of thermal spectra: they are con- 
volutions of convolutions of thermal spectra.l^^lEoSI j^^ Fig. [21] we present the photon 
number spectrum N{E) time-integrated over the 20 s of the entire duration of the 
prompt event of GRB 031203 observed by INTEGRALEU we obtain a typical 
nonthermal power-law spectrum which turns out to be in good agreement with the 
INTEGRAL dataPs^ill and gives evidence of the possibility that the observed GRB 
spectra originate from a thermal emission.'i^ 

Before ending this discussion, we mention that, using the diagrams given in 
Figs. I16H191 in principle one can compute the two free parameters of the fireshell 
model, namely i?*°* and B, from the ratio between the total energies of the P-GRB 
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Fig. 21. Three theo retically predicted time-integrated photon number spectra N(E), computed 
for GRB 031203C2l]are shown here for < t^ < 5 s, 5 < t^ < 10 s and 10 < t;^ < 20 s (dashed and 
dotted curves), where tjj is the photon arrival time at the detector l^^SliMJ The hard to soft behavior 
is confirmed. Moreover, the theoretically predicted time-integrated photon number spectrum N{E) 
corresponding to the first 20 s of the "prompt emission" (black bold curve) is compared with the 
data observed by INTEGRAL.-^ - This curve is obtained as a convolution of 108 instantaneous 
spectra, which are enough to get a good agreement with the observed data. See details in Ref. ll06l 

and of the extended afterglow and from the temporal separation between the peaks 
of the corresponding bolometric light curves. Correspondingly, it is also possible to 
evaluate the temperature at decoupling boosted by the Lorentz gamma factor as well 
as the temperature in the comoving frame of the fireshell. None of these quantities 
depends on the cosmological model. Therefore, one can in principle use this method 
to compute the GRB intrinsic luminosity and temperature at transparency, and 
make GRBs the best cosmological distance indicators. The increase in the number 
of observed sources, as well as the more accurate knowledge of their CBM density 
profiles, will possibly make this procedure viable to test cosmological parameters, 
in addition to the Amati relation.EiMlZl 



4. Progress in theoretical physics 

After describing this general picture of the canonical GRB scenario, it is appropriate 
to turn to some progress in theoretical physics motivated by the quest of clarifying 
some of the main features of the Fireshell model. We will then return to discuss our 
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understanding of some specific GRBs and introduce a new class of GRBs. 

4.1. Progress in understanding the thermalization process of the 
electron-positron-baryon plasma 

I already mentioned in Sec. [T] that a crucial hypothesis underlying the concept of 
the dyadosphere, as well as the dynamical analysis of the PEM and PEMB pulse of 
the fireshell, was the assumption of thermal equilibrium with photons and baryons 
in the e"'"e~ plasma. I also expressed how this assumption has not been adopted 
in alternative models, e.g. the one by Cavallo and Reesp^ where the e+e^ plasma 
degrades due to bremsstrahlung down to temperatures below kT ~ nieC^ and the 
e+e" disappear. The two different concepts of fireball and fireshell then follow. It 
was then necessary to address a theoretical analysis to justify our assumption based 
on first principles. 

In order to establish which scenario is correct we studied the relaxation of an 
electron-positron plasma to thermal equilibrium in Ref. 11181119 1 There the relativis- 
tic Boltzmann equations with exact QED coUisional integrals taking into account 
all relevant two-particle (Compton scattering, Bhabha and Moller scattering, pair 
creation and annihilation in two photons, electron-proton Coulomb scattering) and 
three-particle interactions (relativistic bremsstrahlung, double Compton scattering, 
three photon annihilation and radiative pair creation), see Table [TJ were solved nu- 
merically. It was confirmed that a metastable state called "kinetic equilibrium''^^ 
exists in such a plasma, which is characterized by the same temperature of all parti- 
cles, but nonzero chemical potentials. Such a state occurs when the detailed balance 
of all two-particle reactions is established. It was pointed out in Refs. lll8|119l that 
direct and inverse three-particle interactions are essential in bringing the electron- 
positron plasma to thermal equilibrium. In particular, when inverse three-particle 
interactions are switched off artificially the thermal equilibrium is never reached, 
see Fig. El 

In Ref. 11211 with Aksenov and Vereshchagin we computed numerically the 
relaxation time scales for optically thick electron-positron plasma over a wide 
range of temperatures and baryon loadings using the kinetic code we developed in 
Ref. I118I119I These time scales were previously estimated in the literature by order 
of magnitude arguments using the reaction rates of the dominant processes ! ^ 22x12 3] 
We have shown that these numerically obtained time scales differ from previous 
estimations by several orders of magnitude. 

Temperatures in the range considered in Ref. I118I119I121I 

0.1 < -^ < 10 (12) 

were selected in order to avoid production of other particles such as muons and 
neutrinos. 021 

In the uniform isotropic pair plasma relativistic Boltzmann equation for the 
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Fig. 22. Dependence on time of concentrations of pairs (black), photons (red) and both (thick) 
when all interactions take place (solid). The upper (lower) figure corresponds to the case where 
initially there are mainly pairs (photons). Dotted curves in the upper figure show concentrations 
when inverse triple interactions are neglected. In this case an enhancement of the pairs occurs with 
the corresponding increase in photon number and thermal equilibrium is never reached. 



distribution function fi of the particle species i has the foUowing form: 



1 d 



-iT/.(p.,i) = E('?i-^?/'(p-*))' 



c dt 



(13) 
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Fig. 23. Left: The thermalization time scale of tiie electron-positron-photon component of plasma 
as a function of the total energy density and the baryonic loading parameter. The energy density 
is measured in erg ■ cm~^, time is seconds. Right: The final thermalization time scale of a pair 
plasma with baryonic loading as a function of the total energy density and the baryonic loading 
parameter. The energy density is measured in erg ■ cm~^, time is seconds. 

where the sum is taken over ah two- and three-particle reactions q, and "qf and xi 
are, respectively, the emission and absorption coefficients. 

Table 1. Physical processes included in simulations 



Binary interactions 


Radiative and 

pair producing variants 


M0ller and Bhaba 


Bremsstrahlung 

e+e+ e+ e+ 7 

T + t' +' 

e+ef ^e+ ef -y 


Single Compton 

7eT — >■ 7'eT 


Double Compton 

7eT <-!. -y'e^'-y' 


Pair production 
and annihilation 

77' •<->■ e^e* 


Radiative pair production 
and 3-photon annihilation 

77' <-> 6^6=1=7" 

eTe± <^ 77'7" 


e^7 <->■ e"!" e^ e 


Coulomb scattering 

P1P2 -^ p'iP'2 
peT p'eT 


peT <-^ p'eT 7 
P7 <-4 p'e^e^ 



We developed a relativistic kinetic code for treating the electron-positron-proton 
plasma in the framework of the kinetic Boltzmann equations in the homogeneous 
and isotropic cases.'^^^'^^ We focus only on the electromagnetic interactions. The 
two basic parameters of the plasma are the total energy density p and the baryonic 
loading parameter B = — +''''+ ' where pb, Pj, Pe- i Pe+ are the energy densi- 
ties of baryons, photons, electron, and positrons. We choose the range of plasma 
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parameters relevant to GRB sources 

10^^ <p< lO^^erg • cm-=^ , 10"^ <B<10^ . (14) 

The corresponding temperatures in the thermal equilibrium state are 0.1 < T < 10 
MeV. Firstly, if we neglect neutrino channels as well as the creation and annihilation 
of baryons and the weak interactions, we have no additional particles except for the 
7, e+e~, p. Secondly, given the smallness of the plasma parameter g — (rieA^,)^^, 
where \u is the Debye length, we can use only one-particle distribution functions 
/(p,i) in this most interesting temperature range. The range of the parameter B 
includes cases of almost purely electron-positron plasma, and almost purely electron- 
proton plasma. 

We started from an arbitrary non-equilibrium state for fixed sets of parameters 
p and B. We then solved Eqs. ([T3)) numerically until the plasma steady state is 
reached. As a result, we evaluated the corresponding time scales for the plasma 
relaxation to the equilibrium, reproduced in Fig. [23l We have different times scales 
for the electron-positron component and a time scale for the entire plasma due to 
the large differences in masses of electrons and protons. The final time scales can 
be estimated as rth — inax(T3p,min(Tep,ipp)), where T3p is the three-particle inter- 
action time scale, r^^p is the electron-proton elastic scattering time scale, and tpp 
is the proton-proton elastic scattering time scale, respectively. The exact numeri- 
cal coefficient in the relation between the calculated time scale and the estimated 
formula varies over a wide range of values .'^^B 

The validity of our assumption in the canonical GRB scenario has thus been 
proved. This work can be used as an example of how rigorous fundamental results 
in relativistic plasma physics have been obtained motivated by the study of GRBs. 

4.2. Progress in the physics of neutron stars: Local vs. global 
charge neutrality and critical fields 

Recent progress has been made in the generalization of the Thomas-Fermi model to 
special relativistic regimes as well as to general relativistic regimes (see e.g. Refs.[32l 
I124H127| . We have generalized classical results obtained by Feynman, Metropolis 
and Teller^^Sl and, by the introduction of scaling laws, the classical results obtained 
by Popov and collaborators (e.g. Ref. I129|130p in heavy nuclei to massive cores 
of ~ Mq. On the other hand, we have found special justification when applied to 
astrophysical systems leading to a consistent treatment of white dwarfs (see e.g. 
Ref. l32]) and to a deeper understanding of neutron star physics (see e.g. Ref. 11261 
fT27ll . 

V. S. Popov et al. in Ref. 11241 described heavy nuclei as a degenerate system 
of Nn neutrons, Np protons and N^ electrons constrained to a constant nuclear 
density distribution for the protons and then solved the corresponding relativistic 
Thomas- Fermi equation. We have first generalized the work of V. S. Popov^^^HSS] 
and W. Greiner^^SEMl ^^y eliminating their constraint between the total number of 
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protons and the total number of baryons, Np « A/2, clearly not valid for heavy 
nuclei. We have self-consistently enforced the condition of beta equilibrium in a new 
relativistic Thomas-Fermi equation. Using then the existence of scaling laws these 
results have been extended from heavy nuclei to the case of nuclear matter cores of 
stellar dimensions. In both these treatments a zero Fermi energy of the electrons, 
E^ = 0, was assumed. 

We have recently generalized this dual approach in M. Rotondo et aLli^SI by 
considering first the case of compressed atoms and then, using the existence of scal- 
ing laws, the case of compressed nuclear matter cores of stellar dimensions with a 
positive value of their electron Fermi energies. This approach allows a precise treat- 
ment of the electrodynamical interactions within a compressed atom with all their 
relativistic corrections. Hence, a self-consistent equation of state for compressed nu- 
clear matter is derived which at high densities validates the equation of state due 
to E. E. Salpeter^'^^ and overcomes some of its difficulties, like the appearance of 
negative pressure at low densities. 

We have recently applied such an equation of state to the study of the general 
relativistic white-dwarf equilibrium configurations in J. A. Rueda et al.'^The con- 
tributions of quantum statistics and of the weak and electromagnetic interactions 
have been further generalized there by considering the contribution of the general 
relativistic equilibrium of white dwarf matter. This is expressed by the simple for- 
mula yjgool/^ws — constant, linking the chemical potential of the Wigner-Seitz cell 
^ws with the general relativistic gravitational potential goo at each point of the 
configuration. The configuration outside each Wigner-Seitz cell is strictly neutral 
and therefore no global electric field is necessary to ensure the equilibrium of the 
white dwarf. These equations modify those used by Chandrasekhar by properly ac- 
counting for the Coulomb interaction between the nuclei and the electrons as well as 
inverse beta decay. They also generalize the work of Salpeter by considering a unified 
self-consistent approach to the Coulomb interaction in each Wigner-Seitz cell. The 
consequences for the numerical value of the Stoner-Chandrasekhar-Landau mass 
limit as well as for the mass-radius relation of white dwarfs were then derived.!^ 
This leads to the possibility of a direct confrontation of these results with obser- 
vations. This is currently of great interest because of the cosmological implications 
of the type la supernova^i^SHlMl and because of the low mass white dwarf com- 
panion of the pulsar PSRJ1141-654^i^ as well as the role of white dwarfs in SNe, 
soft gamma-ray repeaters (SGRs) and anomalous X-ray pulsars (AXPs)^^ as an 
explicit alternative to magnetars.l^^ 

These results have been extrapolated to the case of nuclear matter cores of stellar 
dimensions for A « (mpianck/'Tin)'^ ^ 10^^ or Mcore ~ Af© (see Ref. ll25l for details). 
The possibility of obtaining for these systems a self-consistent solution characterized 
by global but not local charge neutrality was explored there. The results generalize 
the considerations presented by Popov et al. in Ref. 11241 corresponding to a nuclear 
matter core of stellar dimensions with zero Fermi energy of the electrons. An en- 
tire family of configurations exists with values of the electron Fermi energy ranging 
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Fig. 24. The electric field in units of the critical field for vacuum polarization Ec = rn^c /{eh) is 
plotted as a function of the coordinate ^, for different values of the electron Fermi energy in units 
of the pion rest energy. The solid line corresponds to the case of zero electron Fermi energy. An 
increase in the value of the electron Fermi energy leads to a reduction in the peak value of the 
electric field. 



from zero to a maximum value {E^)max which is reached when the Wigner-Seitz 
ceh coincides with the core radius. The configuration with E^ = {Ef)„iax corre- 
sponds to the configuration with Np = N,, and Up — 71^: for this hmiting value of 
the Fermi energy the system satisfies both global and local charge neutrality and, 
correspondingly, no electrodynamical structure is present in the core. The other 
configurations generally have overcritical electric fields close to their surface. The 
configuration with E^ = has the maximum value of the electric field at the core 
surface, well above the critical value E^ (see Figs. [53 and [M]). All these cores with 
overcritical electric fields are stable against the vacuum polarization process due to 
the Pauli blocking by the degenerate electrons.'^ 

The above problem is theoretically well defined, and represents a necessary step 
in order to approach the more complex problem of a neutron star core and its 
interface with the neutron star crust. 

Neutron stars are composed of two sharply different components: the liquid 
core at nuclear and/or supra-nuclear densities consisting of neutrons, protons and 
electrons and a crust of degenerate electrons in a lattice of nuclei and possibly 
of free neutrons due to neutron drip when this process occurs (see e.g. Ref. I143p . 
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Fig. 25. The electron Coulomb potential energies in units of the pion rest energy in a nuclear 
matter core of stellar dimensions with A ~ 10^^ or Mcore ~ Mq and Re ~ 10^ cm are plotted as 
a function of the dimensionless variable ^ for different values of the electron Fermi energy, also in 
units of the pion rest energy. The solid line corresponds to the case of zero electron Fermi energy. 
Increasing the value of the electron Fermi energy, the electron Coulomb potential energy depth is 
reduced. 



Consequently, the boundary conditions for the electrons at the surface of the neutron 
star core will generally have a positive value of the electron Fermi energy in order to 
take into account the compressional effects of the neutron star crust on the core (see 
e.g. Ref. I126p . The case of zero electron Fermi energy corresponds to the limiting 
case of the absence of the crust. 

These considerations have been generalized in Ref. ll25l bv looking for a violation 
of the local charge neutrality condition over the entire configuration, still keeping 
its overall charge neutrality. This effect cannot occur locally, and requires a global 
description of the equilibrium configuration. 

I illustrated^^ this novel approach with Rotondo, Rueda and Xue by considering 
the simplest, nontrivial self-gravitating system of degenerate neutrons, protons and 
electrons in beta equilibrium in the framework of relativistic quantum statistics 
and the Einstein-Maxwell equations. The impossibility of imposing the condition 
of local charge neutrality on such systems has been proven in complete generality. 
The crucial role of the constancy of the generalized electron Fermi energy was then 
emphasized and consequently the coupled system of the general relativistic Thomas- 
Fermi equations and the Einstein-Maxwell equations was solved. An explicit solution 
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corresponding to a violation of the local charge neutrality condition over the entire 
star, still satisfying global charge neutrality when electromagnetic, weak and general 
relativistic effects are taken into account, is given there. 

The results presented in Ref. 11251 on nuclear matter cores of stellar dimensions 
evidence the possibility of having the existence of critical electromagnetic fields at 
the core surface. Such an analysis has been extended further by considering the 
case of a neutron star. At nuclear and supranuclear densities a core described by 
a self-gravitating system of degenerate neutrons, protons and electrons within the 
framework of relativistic quantum statistics and Einstein-Maxwell equations was 
considered .'^221 At densities lower than the nuclear density such a core is surrounded 
by a crust. A globally neutral neutron star configuration was examined there in 
contrast with the traditional ones constructed by imposing local charge neutrality. 
To illustrate the application of this approach the Baym, Bethe and PethiclP^ strong 
interaction model of the baryonic matter in the core and in the white dwarf-like 
material of the crust was adopted there. The existence of an overcritical electric 
field at the boundary of the core as predicted in Ref. 11251 was confirmed there. The 
electric field extends over a thin shell of thickness ^ h/{m^c) between the core and 
the crust and becomes largely overcritical in the limit of decreasing values of the 
crust mass (see Refs. 11261 for details). 

All the new gravito-electrodynamical effects discussed here deserve further anal- 
ysis in view of the recent developments in high-energy astrophysics pointing to the 
relevance of overcritical electric fields in neutron stars and black holes .1^ 

The extension of all the above considerations by describing the strong inter- 
action between nucleons through sigma-omega-rho meson exchange in the context 
of the extended Walecka model, all properly expressed within general relativity, 
is currently being developed. In this case it is shown that, exactly as in the non- 
interacting case, the thermodynamic equilibrium condition given by the constancy 
of the Fermi energy of each particle-species can be properly generalized to include 
the contribution of all fields.E^H 

We were clearly motivated in this theoretical development by the study of GRBs. 
The gravitational collapse of a neutron star endowed with a critical electric field 
clearly leads to a violent pair creation process. The gravitational energy extraction 
process due to the presence of an electric field as envisaged in Fig. [Til now becomes 
naturally implementable. The gravitational collapse process is triggered by the neu- 
trons and the protons reaching of a fully relativistic regime. Consequently, their 
gravitational collapse leads to an increase of the charge-to-mass ratio of the core 
and an increase of the electric field given in Fig. [24] well above the critical value 
and over the Pauli blocking of the degenerate electrons of the core. This process is 
currently being studied. Due to the increase of density in the gravitational collapse 
in presence of an overcritical field, some striking analogies have surfaced with the 
physics of ultra high energy heavy ion collisions. 
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4.3. Progress in the physics of dy ado sphere: The dyado-torus of a 
rotating black hole 

In parallel with the above investigations, we have generalized the concept of the 
dyadosphere to the more general Kerr-Newman geometry where the corresponding 
dyadosphere region is instead only axially symmetric. This approach leads to the 
concept of the "dyadotorus," recently introduced in Ref. |4TJ The details have been 
presented in Ref. 11451 The goal was to identify the region in the Kerr-Newman 
geometry where vacuum polarization processes may occur, leading to the creation 
of e~ — e"*" pairs, generalizing the original concept of the "dyadosphere" initially in- 
troduced for Reissner- Nordstrom geometries.^^ ^-^The topology of the axially sym- 
metric dyadotorus was studied for selected values of the electric field, and its elec- 
tromagnetic energy was estimated using three different methods all of which give 
the same result. 

Vacuum polarization processes can occur in the overcritical field of a Kerr- 
Newman black hole inside the dyadotorus. Such a region has an invariant character, 
i.e. its existence does not depend on the observer measuring the electromagnetic 
field: therefore, it is a true physical region. 

Some pictorial representations of the boundary surface similar to those com- 
monly used in the literature have been presented in Ref. 11451 employing Cartesian- 
like coordinates (i.e. interpreting the Boyer-Lindquist radial and angular coordinates 
as ordinary spherical coordinates on flat space) as well as Kerr-Schild coordinates. 
The dyadotorus has been also shown on the corresponding embedding diagram, 
which gives the correct geometry allowing the visualization of spacetime curvature. 

We have then estimated the electromagnetic energy contained in the dyadotorus 
by using three different approaches, all of which give rise to the same final expres- 
sion for the energy. The first one follows the standard approach consisting of using 
the (unnormalized) timelike Killing vector through the Boyer-Lindquist constant 
time slice of the Kerr-Newman spacetime (see e.g. Ref. [82]) . the second one fol- 
lows a recent observer dependent definition by Katz, Lynden-Bell and Bicalii^SUiT] 
for axially symmetric asymptotically flat spacetimes, for which we have used the 
Painleve-Gullstrand geodesic family of infalling observers through the Painleve- 
GuUstrand constant time slice, and the last one adopts the pseudotensor theory (see 
e.g. Ref. I148p . We have found by rough estimates that the extreme Kerr-Newman 
black hole leads to larger values of the electromagnetic energy as compared to a 
Reissner-Nordstrom black hole with the same total mass and charge. 

It is appropriate to recall that the release of energy via the electron-positron pairs 
in the dyadotorus is the most powerful way to extract energy from black holes and 
in every way corresponds to a new form of energy: the "blackholic" energy.!^ This 
is a new form of energy different from the traditional ones known in astrophysics. 
The thermonuclear energy has been recognized to be energy source of main sequence 
stars lasting for 10^ years,'I^the gravitational energy released by accretion processes 
in neutron stars and black holes has explained the energy observed in binary X-ray 
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Fig. 26. The projections of the dyadotorus on the X-Z plane corresponding to different values 
of the ratio \Ei\/Ec = k are shown in the upper panel for /i = 10 and A = 1.49 X 10~^. The 
corresponding plot for the dyadosphere with the same mass energy and charge-to-mass ratio is 
shown in the lower panel for comparison. See details in Ref. 11451 



sources on time scales of 10^ — 10* yearsl^^ The "blackholic" energy appears to be 
the energy source for the most transient and most energetic events in the universe, 
the GRBsP 

We now turn to the understanding of specific GRBs. 
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5. On a new class of "disguised" short GRBs within the Fireshell 
model 

I proceed now to some recent applications of the Fireshell model to infer properties 
of GRBs and of their progenitors. 

In the current literature, "long" GRBs are traditionally related to the idea of a 
single progenitor, identified as a "coUapsar" .^ Similarly, short GRBs are assumed 
to originate from binary mergers formed by white dwarfs, neutron stars, and black 
holes in all possible combinations. It also has been suggested that short and long 
GRBs originate from different galaxy types. In particular, short GRBs are proposed 
to be associated with galaxies with low specific star forming rate.^^^ Some evidence 
against such a scenario have been advanced, due to the small sample size and the 
different estimates of the star forming rates.'i^ 

I give some specific examples of how the understanding of GRB structure and 
of its relation to the CBM distribution, within the fireshell model, leads to a more 
complex and interesting perspective than the one in the "majority view" of the 
current literature. 

In the context of the fireshell model, we have considered a new class of GRBs 
pioneered by Norris and BonnelLli^ This class is characterized by an occasional 
softer extended emission after an initial spike-like emission. The softer extended 
emission has a peak luminosity lower than the one of the initial spike-like emission. 
As shown in the prototypical case of GRB 97022»22lwe developed with Maria Grazia 
Bernardini, and then in GRB Oeoei-Pll and in GRB 071227221 we developed with 
Letizia Caito, we have identified the initial spike-like emission with the P-GRB 
and the softer extended emission with the peak of the extended afterglow (see 
Fig. [27]) . That the time- integrated extended afterglow luminosity is much larger 
than the P-GRB one is crucial. This unquestionably identifies these sources as 
canonical GRBs with B > 10^"* (see Figs. fTBlflT)) . The consistent application of the 
fireshell model allows us to infer the CBM filamentary structure and average density, 
which, in that specific case, is ricbm ^ 10^'' particles/cm'^. ^° This low CBM density 
value explains the peculiarity of the low extended afterglow peak luminosity and its 
more protracted time evolution (see Fig. [28]). These features are not intrinsic to the 
progenitor, but depend uniquely on the peculiarly low value of the CBM density. 

This led us to expand the traditional classification of GRBs to three classes: 

• "Genuine" short GRBs: These are the ones in which the total energy emitted 
in the P-GRB is greater than the one emitted in the extended afterglow. They 
are characterized by a small value of the fireshell baryon loading B < 10~^. 

• "Disguised" short GRBs: These are the ones in which the total energy 
emitted in the P-GRB is smaller than that emitted in the extended afterglow, 
but the P-GRB peak luminosity is greater than that of the extended afterglow. 
They are characterized by a high value of the fireshell baryon loading B > 3.0 x 
10~^, and a low value of the CBM average density ricbm ~ 10"'^ particles/cm^ 
which "deflates" the extended afterglow peak luminosity. 
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Fig. 27. The "canonical GRB" light curve computed theoretically for the prompt emission of 
GRB 970228. The BeppoSAX GRBM (40-700 keV, above) and WFC (2-26 keV, below) light 
curves (data points) are compared with the extended afterglow peak theoretical ones (solid red 
lines). The onset of the extended afterglow coincides with the end of the P-GRB (represented 
qualitatively by the dotted blue lines). For this source we have B ~ 5.0 x 10"'^ and {n^i,^) ~ 10~^ 
particles/cm^. The total time-integrated energy emitted in the P-GRB is therefore much lower than 
the one emitted in the extended afterglow, but this last one has its peak luminosity "deflated" by 
the low CBM density. See details in Ref. |90l 



• The remaining "long duration" ones: These are the ones in which the 
total energy emitted in tlie P-GRB is smaller than that emitted in the extended 
afterglow, and the P-GRB peak luminosity is smaller than that of the extended 
afterglow. They are characterized by a high value of the fireshell baryon loading 
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Fig. 28. The theoretical fit of the BeppoSAX GRBM observations (solid red line, see Fig. [27ll 
is compared with the extended afterglow light curve in the 40-700 keV energy band obtained 
rescaling the CBM density to (ricim) = 1 particle/cm^ keeping constant its shape and the values 
of the fundamental parameters of the theory B*°' and B (black double-dotted line). The P-GRB 
duration and luminosity (blue dotted line), depending only on _E*°' and B, are not affected by 
this process of rescaling the CBM density. See details in Ref. 1901 



B > 3.0 X 10 ■*, and a high value ol the CBM average density nchm 
particle/cm'^. 



> 1 



A CBM density Ucbm ^ 10^^ particles/cm'^ is typical of a galactic halo environ- 
ment, and GRB 970228 -was indeed found to be in the halo of its host galaxy.E^^li^ 
We therefore proposed that the progenitors of this ne-w class of "disguised" short 
GRBs are merging binary systems, formed by neutron stars and/or ■white dwarfs 
in all possible combinations, which spiraled out from their birth place into the 
jjg^jQJ90lT53lT57| This hypothesis can also be supported by other observations (e.g. in 
the optical band), clearly showing an offset of the GRB position from the center 
of the host galaxy (see e.g. the case of GRB 050509tP^. Vice versa, there is also 
evidence for the opposite correlation, directly relating the soft tail peak luminosity 
to the CBM density: GRBs displaying a more luminous prolonged soft tail appear 
to have a systematically smaller offset from the center of their host galaxy .I^SMlfiU 

I turn now to some specific examples of "disguised" short GRBs. 



6. Application to GRB 970228 

The GRB 970228 has been the first "afterglow" ever detected.'^ It was quite surpris- 
ing that, as we returned to its analysis years later with M.G. Bernardini in Ref.f90l 
it was indeed discovered that this GRB became the prototype for the disguised 
short GRB class identified by Norris & Bonnellfi^ 
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The GRB 970228 was detected by the Gamma-Ray Bm^st Monitor (GRBM, 
40-700 keV) and Wide Field Cameras (WFC, 2-26 keV) on board BeppoSAX on 
February 28.123620 UTP^The burst prompt emission is characterized by an initial 
5 s strong pulse followed, after 30 s, by a set of three additional pulses of decreasing 
intensity P^ Eight hours after the initial detection, the NFIs on board BeppoSAX 
were pointed at the burst location for a first target of opportunity observation and 
a new X-ray source was detected in the GRB error boxPil A fading optical transient 
has been identified in a position consistent with the X-ray transient ,1^^ coincident 
with a faint galaxy with redshift z ~ 0.695. ^^'^ Further observations by the Hubble 
Space Telescope clearly showed that the optical counterpart was located in the 
outskirts of a late-type galaxy with an irregular morphology .1^^ 

The BeppoSAX observations of the GRB 970228 prompt emission revealed a 
discontinuity in the spectral index between the end of the first pulse and the begin- 
ning of the three additional ox\b^^^^^^ The spectrum during the first 3 s of the 
second pulse is significantly harder than during the last part of the first pulse,'^^2nS2l 
while the spectrum of the last three pulses appear to be consistent with the late 
X-ray afterglow.'^^^'^^ This was soon recognized by Ref. I110|162l as pointing to an 
emission mechanism producing the X-ray afterglow already taking place after the 
first pulse. 

As recalled above, the simultaneous occurrence of an extended afterglow with 
total time-integrated luminosity larger than that of the P-GRB, but with a smaller 
peak luminosity, was indeed explainable in terms of a peculiarly small average value 
of the CBM density and not due to the intrinsic nature of the source. We have 
showrpSlthat GRB 970228 is a very clear example of this situation. We identify the 
initial spikelike emission with the P-GRB, and the late soft bump with the peak of 
the extended afterglow. GRB 970228 shares the same morphology and observational 
features with the sources analyzed by Norris & Bonnell,'^^ e.g. GRB 050709,*^^ 
GRB 050724i£Slas weh as GRB 060614.ESi Therefore, we proposed GRB 970228 as 
the prototype for this new GRB class.l^ 



6.1. The analysis of the GRB 970228 prompt emission 

In Fig. [27] I recall the theoretical fit of the BeppoSAX GRBM (40-700 keV) and 
WFC (2-26 keV) light curves of the GRB 970228 prompt emissiorii^^ performed by 
Bernardini et al.!^ We identified the first main pulse with the P-GRB and the three 
additional pulses with the extended afterglow peak emission, consistent with the 
above mentioned observations by Ref. [51] and Ref. 11621 The last three such pulses 
have been reproduced assuming three overdense spherical CBM regions (see Fig. [29ll 
with a very good agreement (see Fig. [27| . 

We therefore obtain for the two parameters characterizing the source in our 
model Elf = 1.45 x lO^'' erg and B = 5.0 x lO-^l^n! This implies an initial e+e" 
plasma created between the radii ri — 3.52 x 10^ cm and r2 = 4.87 x 10® cm with 
a total number of e+e^ pairs Ng± = 1.6 x 10^^ and an initial temperature T = 1.7 
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Fig. 29. The CBM density profile we assumed to reproduce the last three pulses of the 
GRB 970228 prompt emission (red line), together with its average value {n^i,jy^) = 9.5 X 10~* 
particles/cm"^ (green line). Details in Ref. 1901 

jyjgYJM! The theoretically estimated total isotropic energy emitted in the P-GRB 
is Ep-oRB — l.l%i?*±* = 1.54 X 10^^ erg, in excellent agreement with the energy 
observed in the first main pulse (-Ep^cfls ^ ^-^ ^ ^^^^ ^^S in 2 — 700 keV energy 
band, see Fig. [27l), as expected due to their identification.ISS 

After the transparency point at ro = 4.37 x 10^'* cm from the progenitor, the 
initial Lorentz gamma factor of the fireshell is 70 = 199. On average, during the 
extended afterglow peak emission phase we have for the CBM (TZ) = 1.5 x 10~^ and 
{ncbm) — 9.5 X 10~^ particles/cm^. This very low average value for the CBM density 
is compatible with the observed occurrence of GRB 970228 in its host galaxy's 
jjg^jJi55lT56|I62l a,nd it is crucial in explaining the light curve behavior. 

The values of -B*±' and B we determined are univocally fixed by two tight con- 
straints. The first one is the total energy emitted by the source all the way up to 
the latest afterglow phases (i.e. up to ~ 10® s). The second one is the ratio between 
the total time-integrated luminosity of the P-GRB and the corresponding one of 
the whole extended afterglow (i.e. up to ^ 10® s). In particular, in GRB 970228 
such a ratio turns out to be ^ 1.1%, typical of a "long" GRB (see Figs. [T6lfT7| . 



6.2. Rescaling the CBM density 

I recall now an explicit example in order to probe the crucial role of the average 
CBM density in explaining the relative intensities of the P-GRB and of the extended 
afterglow peak in GRB 970228 (for details see Ref. |90]). We keep fixed the basic 
parameters of the source, namely the total energy E"*"' and the baryon loading B, 
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therefore keeping fixed the P-GRB and the extended afterglow total time- integrated 
luminosities. Then in a gedanken experiment we rescale the CBM density profile 
given in Fig. [511 by a constant numerical factor in order to raise its average value 
to the standard one {ricbm) = 1 particle/cm^. We then compute the corresponding 
light curve, shown in Fig. [28l 

We notice a clear enhancement of the extended afterglow peak luminosity with 
respect to the P-GRB one, when compared to and contrasted with the observational 
data presented in Fig.[27l The two light curves actually cross at if ~ 1.8 x 10'* s since 
their total time-integrated luminosities must be the same. The GRB "rescaled" to 
(ncbm) = 1 particle/cm'^ appears to be totally similar to, e.g., GRB 050315^^ and 
GRB 991216.I S3I93 I 101 ] 

It is appropriate to emphasize that, although the two underlying CBM density 
profiles differ by a constant numerical factor, the two extended afterglow light curves 
in Fig. [28] do not. This is because the absolute value of the CBM density at each 
point affects in a nonlinear way all the subsequent evolution of the fireshell due to 
the feedback on its dynamics.'^ Moreover, the shape of the surfaces of equal arrival 
time of the photons at the detector (EQTS) is strongly elongated along the line of 
sight.'^S Therefore, this is a good example of how photons coming from the same 
CBM density region are observed over a very long arrival time interval during the 
extended afterglow. 



6.3. GRB 970228 and the Amati relation 

It is appropriate now to turn to the "Amati relation"ESSIiZ2] between the isotropic 
equivalent energy emitted in the prompt emission Eiso and the peak energy of 
the corresponding time-integrated spectrum Ep^i in the source rest frame. It has 
been shown by Ref. I169I170I that this correlation holds for almost all the "long" 
GRBs which have a redshift and a measured Ep^i , but not for the ones classified 
as "short" .11221 

It clearly follows from our treatment that for the "canonical GRBs" with large 
values of the baryon loading and high {ricbm), which presumably are most of the 
GRBs for which the correlation holds, the leading contribution to the prompt emis- 
sion is the extended afterglow peak emission. The case of the "fake" short GRBs is 
completely different: it is crucial to consider both components separately since the 
P-GRB contribution to the prompt emission in this case is significant. 

To test this scenario, we evaluated from our fit of GRB 970228 Eiso and Ep^i 
only for the extended afterglow peak emission component, i.e. from tji^ = 37 s to 
f^ = 81.6 s. We found an isotropic energy emitted in the 2-400 keV energy band 
E^so = 1-5 X 10^2 erg, and Ep^, = 90.3 keV^ As is clearly shown in Fig. [30l the 
sole extended afterglow component of GRB 970228 prompt emission is in perfect 
agreement with the Amati relation.E^ 

From the energetics of this source, as well as from its location in the galactic 
halo and the low density of the CBM, we conclude that the progenitors of GRB 
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Fig. 30. The estimated values for Ep^i and Eigo obtained by our analysis (black dot) compared 
with the "Amati relation". 1^^ the solid line is the best fitting a power la W^ ^^ la nd the dashed lines 
delimit the region corresponding to a vertical logarithmic deviation of 0.4l ^^^ lThe uncertainty in 
the theoretical estimated value for Epi has been assumed conservatively to be 20%. See details in 

Ref-irm 



970228 are merging binary neutron stars. 



7. Application to GRB 060614 

I turn now to GRB 0606l4 ^^ l ^^^l for at least three different reasons. 

The first novehy is that it is the first clear example of a nearby (2 = 0.125) 
long GRB not associated with a bright Ib/c supernova (SN).E^^E21! Jt has been 
estimated that, if present, the SN-component should be about 200 times fainter 
than the archetypal SN 1998bw associated with GRB 980425; moreover, it would 
also be fainter (at least 30 times) than any stripped-envelope SN ever observed .^^^ 

Within the standard scenario, long duration GRBs (Tgo > 2 s) are thought 
to be produced by SN events during the collapse of massive stars in star forming 
regions ("coUapsar").^- The observations of broad- lined and bright type Ib/c SNe 
associated with GRBs are often reported to favor this scenario (see Ref. 77. and 
references therein). The ansatz of the coUapsar model has been that every long 
GRB should have a SN associated with it LUSI Consequently, in all nearby long GRBs 
(2: < 1), SN emission should be observed. 

For these reasons the case of GRB 060614 is of great relevance. Some obvious 
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hypotheses have been proposed and ruled out: the chance superposition with a 
galaxy at low redshifiji^ and strong dust obscuration and extinctionP^ Appeal 
has been made to the possible occurrence of an unusually low luminosity stripped- 
envelope core-collapse SN.^^'^ 

The second novelty of GRB 060614 is that it challenges the traditional separation 
between long soft GRBs and short hard GRBs. TraditionallyJi^SEini the "short" 
GRBs have Tgg < 2 s, present an harder spectrum and negligible spectral lag, 
and are assumed to originate from the merging of two compact objects, i.e. two 
neutron stars or a neutron star and a black hole (see e.g. Refs. l60l68l87ll80Hl82l and 
references therein). GRB 060614 lasts about one hundred seconds (Tgo = (102 ± 5) 
g^ 1166] it fulfills the Ep^^^-Eiso correlation,'^^ and therefore traditionally it should be 
classified as a "long" GRB. However, its morphology is different from typical long 
Qj^gg[T76liMl Jts optical afterglow luminosity is intermediate between the traditional 
long and short ones.'^^ Its host galaxy has a moderate specific star formation rate 
(RHost « 2M,y-\L*)-\ M^Host « -15.5).E2ll23 The spectral lag in its light 
curves is very small or absent.'^^ All these features are typical of short GRBs. 

A third novelty of GRB 060614 is that its 15-150 keV light curve presents a 
short, hard and multi-peaked episode (about 5 s). The episode is followed by a softer, 
prolonged emission that manifests a strong hard to soft evolution in the first 400 s of 
data.li^IlThe total fluencc in the 15-150 keV energy band is F = (2.17±0.04) x 10"^ 
erg/cm^, the 20% emitted during the initial spikelike emission, where the peak 
luminosity reaches the value of 300 keV before decreasing to 8 keV during the 
BAT-XRT overlap time (about 80 s). 

These apparent contradictions find a natural explanation in the framework of the 
"fireshell" model, as explicitly shown with L. Caito.li^ Within the fireshell model, 
the occurrence of a GRB-SN is not a necessity. 



7.1. The fit of the observed luminosity 

We have proceeded to interpret GRB 060614 as a "disguised" short GRB.EsU We 
have performed the analysis of the observed light curves in the 15-150 keV energy 
band, corresponding to the 7-ray emission observed by the BAT instrument on 
the Swift satellite, and in the 0.2-10 keV energy band, corresponding to the X-ray 
component from the XRT instrument on the Swift satellite. The optical emission 
represents less than 10% of the total energy of the GRB and is therefore neglected. 
From this fit (see Fig. [STjl we have derived the total initial energy E^^^ , the value 
of B as well as the effective CBM distribution (see Fig.ES]). We find i^fjj = 2.94 x 
10^^ erg, which accounts for the bolomctric emission of both the P-GRB and the 
extended afterglow. Such a value is compatible with the observed Eiso — 2.5 x 10*^^ 
ej-gliSfil The value oi B is B = 2.8 x 10"'^ corresponds to a canonical GRB with a 
very clear extended afterglow which energetically dominates the P-GRB. From the 
model, having determined E^^^ and B, we can compute the theoretically expected 
P-GRB energetics Ep-cjiB-^e obtain Ep-gb.b — 1-15 x 10^" erg, which is in 
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good agreement with the observed Eiso,ip — 1-18 x 10^° ergJiSSIxhe Lorentz gamma 
factor at the transparency point is 70 = 346, one of the highest of all the GRBs we 
have examined. 

In Fig. [3T] we plot the comparison between the BAT observational data of the 
GRB 0606014 prompt emission in the 15-150 keV energy range and the P-GRB and 
extended afterglow light curves computed within our model. The temporal variabil- 
ity of the extended afterglow peak emission is due to the inhomogeneities in the 
effective CBM density (see Figs. EH [32|). Toward the end of the BAT light curve, 
the good agreement between the observations and the fit is affected by the Lorentz 
gamma factor decrease and the corresponding increase of the maximum viewing 
angle. The source visible area becomes larger than the typical size of the filaments. 
This invalidates the radial approximation we use for the CBM description. To over- 
come this problem it is necessary to introduce a more detailed three-dimensional 
CBM description, in order to avoid an over-estimated area of emission and, cor- 
respondingly, to describe the sharpness of some observed light curves. We are still 
working on this issue. *^ 186 188 

We turn now to the crucial determination of the CBM density, which is de- 
rived from the fit. At the transparency point it has the value ricbm = 4.8 x 10^^ 
particles/cm"^ (see Fig. [5^ . This density is compatible with the typical values of the 
galactic halos. During the peak of the extended afterglow emission the effective av- 
erage CBM density decreases reaching (ricbm) ~ 2.25 x 10^^ particlcs/cm"^, possibly 
due to an ongoing fragmentation of the shelP^ or due to a fractal structure in the 
CBM. The TZ value on average was (TZ) = 1.72 x 10~^. Note the striking analogy of 
the numerical value and the overall radial dependence of the CBM density in the 
present case of GRB 060614 when compared and contrasted with the ones of GRB 
970228.'S2l More details can be found in Ref. IT531 

In view of the low CBM density and the overall energetics, which is considerably 
smaller than for GRB 970228, we conclude that the progenitor of GRB 060614 is a 
binary merger formed by a neutron star and a white dwarf. Both GRB 970228 and 
GRB 060614 are clear counterexamples to the proposal that binary mergers lead 
necessarily to short GRBs.'^^M] j^^ both cases they lead to long GRBs disguised as 
short ones. 



8. GRB 071227: an additional case of a disguised short burst 

GRB 071227 presents some intriguing anomalies and has been shown to represent 
a third case of a disguised short GRB.HSU As for the "Norris and Bonncl" GRBs, 
its BAT light curve shows in the 15-150 keV range a multi-peaked structure lasting 
Tgo = (1.8 ± 0.4) s, followed by an extended but much softer emission up to to + 100 
g[l89] ^ fading X-ray (0.3-10 keV) and a faint optical afterglow have also been 
identified. The optical afterglow emission allowed the measurement of its redshift, 
z — 0.383, and therefore of its isotropic equivalent energy, Eiso — 5.8 x 10^° erg in 
20-1300 keV.'^^ The observed X-ray and optical afterglow is superimposed on the 
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Fig. 31. The BAT 15—150 keV light curve (green points) at 1 s time resolution compared with the 
corresponding theoretical extended afterglow light curve we compute (red line). The onset of the 
extended afterglow is at the end of the P-GRB (qualitatively sketched in blue lines and delimited 
by dashed gray vertical lines). Therefore the zero of the temporal axis is shifted by 5.5 s with 
respect to the BAT trigger time. The peaks of the extended afterglow light curves are labeled to 
match them with the corresponding CBM density peak in Fig. 1321 In the upper right corner there 
is an enlargement of the P-GRB at 50ms time resolution (reproduced from Ref. I172II showing its 
structure. See details in Ref. 11531 



plane of the host galaxy, at (15.0 ± 2.2) kpc from its center. 

On the basis of these characteristics, GRB 071227 has been classified as a short 
burst. This statement is supported by other main features: 1) If we consider the first 
and apparently predominant short-duration episode, it does not satisfy the Amati 
relation between the isotropic equivalent radiated energy of the prompt emission 
Eiao and the cosmological rest- frame vF^, spectrum peak energy E^^i?-^^ ^^" ^^■^ ^^'^ 
2) The spectral lag of the first spike-like emission in the 25-50 keV to 100-350 keV 
bands is consistent with zero.'^^ 3) Multiwavelength observations performed over 
many days have shown that there is no association with a Ib/c hypernova, the type 
of SN generally observed with GRBs, even if it is a nearby burst and its isotropic 
energy is compatible with that of other GRBs associated with them,'iS21 although 
the upper limits are not deep enough to rule out a low-energetic core-collapse event. 
Nevertheless, the explosion of this burst in a star-forming region of a spiral galaxy 
and its prolonged tail of emission make it most likely to be a long burst. 

With L. Caito^^ we show that all these ambiguities and peculiarities can be 
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Fig. 32. The effective CBM density (red line) and the TZ parameter (blue hne) versus the radial 
coordinate of the shell. The CBM density peaks are labeled to match them with the corresponding 
extended afterglow light curve peaks in Fig. 1311 They correspond to filaments of characteristic size 
Ar ~ 10"'^^ cm and density contrast A?icbm/("ci)m) ^ 20 particles/cm"^. See details in Ref. 11531 

explained in the framework of the fireshell model if we assume GRB 071227 to be 
a disguised short burst, in which the first spike-like emission coincides with the P- 
GRB and the prolonged softer tail with the peak of the extended afterglow emitted 
in a low CBM density region. We show, moreover, that this tail satisfies the Amati 
relation, and this is consistent with our interpretation. 



8.1. The interpretation of the GRB 071227 light curves 

We have analyzecP^ the observed light curves of this burst in the 15-50 keV band- 
pass, corresponding to the lowest band of the gamma-ray emission, detected by 
the BAT instrument on the Swift satellite, and in the 0.3-10 keV energy band, 
corresponding to the X-ray component from the XRT instrument. To model the 
CBM structure, we assume that ricbm is a function only of the radial coordinate, 
ncbm — ^cbmiT) (radial approximation). The CBM is arranged in spherical shells 
of width ~ IQ-'^^-lO-'^^ cm arranged in such a way that the corresponding modula- 
tion of the emitted flux closely resembles the observed shape. We assumed that the 
first short spike-like emission represents the P-GRB and the gamma-ray tail is the 
peak of the extended afterglow. We therefore began the simulation in such a way 
that the extended afterglow light curve begins in coincidence with the peak of the 
P-GRB (about Is), as shown in Fig. [33l To reproduce the observational data and 
the energetics observed for the P-GRB emission {Eiso ^ 1-0 x 10^"'^ erg), we require 
the initial conditions E^^^ = 5.04 x 10^^ erg and B = 2.0 x 10"''. In Figs. [33] and 
[341 we plot the comparison of the GRB 071227 BAT and XRT data with the the- 
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Fig. 33. The BAT 15-50 keV light curve (green points) compared with the corresponding theo- 
retical extended afterglow light curve (red line). The P-GRB is qualitatively sketched by the blue 
line. See details in Ref. 11541 



oretical extended afterglow light curves. We obtained a good result in the prompt 
emission for the 15-50 keV bandpass (see Fig. [33|) . while for 0.3-10 keV we only suc- 
ceeded in reproducing the first decaying part of the XRT light curve (see Fig. |34|. 
We assumed this to correspond to the possible onset of the "plateau" phase of the 
extended afterglow.^^^ 

From our simulation,'-*^ the amount of energy stored in the P-GRB is found 
to be about 20% of the total energetics of the explosion. Hence, this burst cannot 
be a short burst within the fireshell scenario. The baryon loading obtained {B = 
2.0 X 10"'*) remains in the range of long duration GRBs. This is a very critical value, 
because it is very close to the crossing point of the plot of the energetics of GRBs as a 
function of B (see above). This is the lowest baryon loading that we have ever found 
in our analysis within the fireshell scenario. From our analysis, we found a peculiar 
result for the average CBM density. We obtained a density of ricbm = 1.0 x 10~^ 
particles/cm"^ at the beginning of the process, later decreasing to ricbm — 1-0 x 10~* 
particles/cm"^. This low average density, inferred from the analysis, is responsible for 
the strong deflation of the gamma-ray tail. However, at the radius of about 2.0 x lO^'' 
cm, the density becomes higher and reaches the value of ricbm = 10 particles/cm^ 
(the complete profiles of ricbm and TZ as functions of the radial coordinate are given 
in Fig. [35|) . This is compatible with the observations. The observed X-ray and 
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Fig. 34. The XRT 0.3-10 keV light curve (green points) compared with the corresponding theo- 
retical extended afterglow light curve we obtain (red line). The X-ray data corresponding to the 
first 60 seconds are not available since XRT starts to observe only 60 seconds after the BAT trigger. 
We stop our analysis at 7 ~ 5 when our relativistic dynamical model can no longer be applied. 
See details in Ref. Il54i 



optical afterglow of GRB 071227 is indeed superimposed on the plane of the host 
galaxy, at (15.0 ±2.2) kpc from its center.ESSl An interesting possibility observed by 
D. Arnett (private communication) is that this very low density "cavity" could be 
formed in the coalescing phase of a binary formed by a neutron star and a white 
dwarf. Accurate studies of compact object mergers have shown the distribution 
of merger locations for different host galaxies. ^^^^^ In starburst galaxies, most 
of the mergers are expected to be found within hosts, while in elliptical galaxies 
a substantial fraction of mergers take place outside hosts. Spiral galaxies, hosting 
both young and old stellar populations, represent the intermediate case between the 
preceding two. This result is therefore compatible with our hypothesis about the 
binary nature of the progenitor of GRB 071227. Although they did not consider the 
case of binary systems formed by a neutron star and a white dwarf, the progenitor 
of GRB 071227 would fit into the tight binary scenario described in Ref. 11931 

These results clearly imply that GRB 071227 is another example of a disguised 
short burst.E^ 



July 6, 2011 0:9 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



55 




0.0x10° 5.0x10^^ 1.0x10^^ 1.5x10^^ 2.0x10^^ 
Radial coordinate (cm) 

Fig. 35. The CBM particle number density ricbm (upper panel, red line) and the 7^ parameter 
(lower panel, blue line) as functions of the radial coordinate. See details in Ref. 11541 



8.2. GRB 071227 within the Amati relation 

We studied the position of GRB 071227 in the Ep.i-Eiso plane.^^^ At the observed 
redshift, we assumed a "flat A-CDM model" with H^ = 70 Km/s/Mpc and ^\ = 
0.73. For the first, hard spike, lasting about 1.8s, using Konus/WIND data,n^and 
integrating between 1 and 10 000 keV, we found that Ep^i = (1384 ± 277) keV and 
Eiso — (1.0 ±0.2) X 10^^ erg. As shown in Fig.[36l this is inconsistent with the Ep^i- 
Eiso relation, and instead occupies the short-populated region of the plane. For the 
long tail, lasting about 100s, we used a band model with a = —1.5 and fi — —3, 
which are typical of soft events. These values are compatible with the low quality 
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Fig. 36. Location of the initial short spike and soft long tail of GRB 071227 in the Ep^i-Eigo 
plane. The data points of long GRBs are from Ref. 111611901 the data points and limits of short 
GRBs are from Ref. I17UI190I195I The continuous lines show the best-fit power law and the 2a 
confidence region of the correlation, as determined by Ref. 11161 See details in Ref. 11541 



Statistics of this event. We found that Ep,i = 20ti^ kcV and E^so = (2.2±0.1) x 10^^ 
erg. With these values, the tail of emission is fully consistent with the Amati relation, 
as for any long GRB (see Fig. [351) • This clearly supports our hypothesis about the 
nature of GRB 071227.[i^ 



9. Analysis of the "short" GRB 050509b 

I would like to stress, after the considerations we have presented for the last three 
GRB sources, the consequences for GRB 050509b which was heralded with great 
prominence as a short GRB on the cover page of Nature!^^ This is a particularly 
exciting system, since the previous work on the disguised short GRBs can be used 
to properly identify the nature of this intriguing system. ^^^ 

The ob servations of GRB 050509b by BAT and XRT on board the Swift satel- 
lite^SSESSl have represented a new challenge to the classification of GRBs as long 
and short, since it is the first "short" GRB associated with an afterglow .'^^Sl Its 
prompt emission observed by BAT lasts 40 milliseconds, but it also has an after- 
glow in the X-ray band observed by XRT, which begins 100 seconds after the BAT 
trigger (the time needed to point XRT at the position of the burst) and lasts un- 
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til « 1000 seconds. It is located 40 kpc away from the center of its host galaxy 
(see Fig. [571 and Ref. 158) . which is a luminous, non-star- forming elliptical galaxy 
with redshift z = 0.225.^^^ Although an extensive observational campaign has been 
performed using many different instruments, no convincing optical-infrared can- 
didate afterglow nor any trace of any supernova has been found associated with 
GRB OSOSOOb.E^SlSnnHinil An upper limit in the i?-band 18.5 days after the event 
onset implies that the peak flux of any underlying supernova should have been ~ 3 
mag fainter than the one observed for the type Ib/c supernova SN 1998bw asso- 
ciated with GRB 980425, and 2.3 mag fainter than a typical type la supernova 
(Ref. 12031 see also Ref. [202]). An upper limit to the brightening caused by a su- 
pernova or supernova-like emission has also been established at 8.17 days after the 
GRB: Re ^ 25.0 mag.l^^ While some core-collapse supernovae might be as faint as 
(or fainter than) this limit ,'222' the presence of this supernova in the outskirts of an 
elliptical galaxy would be truly extraordinary.'2532Ml 

Unfortunately, we cannot obtain exhaustive observational constraints for this 
GRB because XRT data are missing in between the first 40 milliseconds and 100 
seconds. However, this makes the theoretical work particularly interesting. With 
G. De Barr 0^^221 .^^g have inferred from first principles some characteristics of the 
missing data which are inferred by our model, and consequently we reached an 
understanding of the source. This is indeed the case, specifically, for the verification 
of the Amati relatiorp22llMlMl for these sources. 

We have proceeded with the identification of the two basic parameters, B and 
ncBM, within two different scenarios. '^^ We first investigate the "ansatz" that this 
GRB is the first example of a "genuine" short burst. After disproving this possibility, 
we show that this GRB is indeed another example of a disguised short burst. 

9.1. Data analysis of GRB 050509b 

9.1.1. Scenario 1 

We first attempted to analyze GRB 050509b under the scenario that assumes it is a 
"genuine" short GRB, namely a GRB in which more than 50% of the total energy 
is emitted in the P-GRB. This would be the first example of an identified "genuine" 
short GRB. 

Within our model, the only consistent solution that does not contradict this 
assumption leads to the interpretation that all the data belongs to the extended 
afterglow phase; the BAT data of the prompt emission (see Fig. 2 in Ref. I196P 
are then the peak of the extended afterglow, and the XRT data represents the 
decaying phase of the extended afterglow (which in the literature is simply called 
"the afterglow," see above. 

In Fig. [311 we recall the result of this analysis. We obtained the following set 
of parameters: Sf* = 2.8 x 10^^ erg, B = 1.0 x 10"*, and ucbm = 1-0 x lO'^ 
particles/cm^. These parameters would imply, however, that the energy emitted in 
the P-GRB should be almost 72% of the total value. This P-GRB should have been 
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Fig. 37. Keck LRIS G-band image, zoomed to show the XRT error circle. The larger, blue circle 
is the revised XRT position from Ref. 12091 the smaller, green circle to the west and north of that 
is the 2a confidence region of the XRT position computed in Ref. 11581 The 11 sources consistent 
with the Ref. 12091 X-ray afterglow localization are labeled in the image. North is up and east is to 
the left. Gl is the large galaxy to the west and south of the XRT. Bad pixel locations are denoted 
with "BP". Figure reproduced from Ref. 1 1581 with the kind permission of J. Bloom. 

clearly observable, and has not been detected. Consequently, this scenario is ruled 
out and we conclude that GRB 050509b cannot be interpreted as a "genuine" short 
GRB. 



9.1.2. Scenario 2 

We then analyzed GRB 050509b under the alternative scenario that assumes the 
energy of the extended afterglow is higher than the P-GRB one. Within our model, 
the only consistent solution leads to the prompt emission observed by BAT (see 
Fig. 2 in Ref. ll96p being interpreted as the P-GRB, and the X-ray decaying afterglow 
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Fig. 38. Our numerical simulation within scenario 1, assuming that GRB 050509b is a "genuine" 
short GRB, i.e. that the P-GRB is energetically predominant over the extended afterglow. The BAT 
data (crosses) are interpreted as the peak of the extended afterglow. In this case, the predicted P- 
GRB (solid rectangle) total energy is more than twice that of the extended afterglow. The solid line 
is the theoretical light curve in the 15-150 keV energy band, and the dashed one is the theoretical 
light curve in the 0.3—10 keV energy band. The dot-dashed horizontal line represents the BAT 
threshold and the solid horizontal one represents the XRT threshold. See details in Ref. 11971 



data observed by XRT being interpreted as the extended afterglow. 

In Fig. [39l we show the result of this analysis. We obtained the parameters 
E^^^ — 5.52 X 10"** erg, B = 6x 10"'*, and an almost constant CBM density ucbm = 
1.0 X lO"'^ particles/cm^. The low value of the number density is justified by the 
GRB being located 40 kpc away from the center of the host galaxy (see Fig. [37]).'^^ 
The P-GRB has an estimated energy of Ep-qrb = 28%i?fQj, which means that 
72% of the energy is released in the extended afterglow. The peak of the extended 
afterglow, theoretically predicted by our model in Fig.[39l was not observed by BAT, 
since the energy was below its threshold, and also not observed by XRT, since 
unfortunately its data collection started only 100 seconds after the BAT trigger. 
More details regarding the fulfillment of the Amati relation by this source can be 
found in Ref. [T97l 

Following our classification, therefore, due to the values of the baryon loading 
and of the CBM density, as well as due to the offset with respect to the host galaxy, 
GRB 050509b is consistent with being another example of a disguised short GRB. 
This follows the previous identification of GRB 970228,1221 GRB OeOGl^Sland GRB 
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Fig. 39. Our numerical simulation within scenario 2, assuming that the extended afterglow is 
energetically predominant over the P-GRB. In this case, the predicted P-GRB (solid rectangle) 
is less then twice the extended afterglow. We interpret the BAT data (crosses) as the P-GRB 
and the XRT data as the extended afterglow. The P-GRB has just 28% of the total energy. The 
double-dashed line is the theoretical light curve in the band 15—150 keV, and the dot-dashed line 
is the theoretical light curve in the band 0.3-10 keV. The two horizontal lines are from above to 
below: the BAT threshold and the XRT threshold. See details in Ref. [T97l 



071227.ES11 



9.1.3. Conclusions on GRB 050509b 
We can then conclude that: 



(1) GRB 050509b is not a genuine short burst, but instead a disguised short one, 
namely a canonical long GRB. 

(2) In view of the very low density of the CBM and of the clear location outside 
the host galaxy as well as from the very stringent observational limits on the 
absence of any supernova, we can safely conclude that the progenitor of this 
system is a binary merger possibly formed by a neutron star and a white dwarf 
or possibly two white dwarf components. 123 

(3) Once again, this is an additional example of a GRB not originating from a 
coUapsar. 
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10. General conclusions for the disguised short GRBs 

We have just shown how the existence of the disguised short GRBs gives a coun- 
terexample both to the coUapsar model and to the proposal that binary mergers 
lead necessarily to short GRBs. It is interesting that the above considerations fol- 
lowing from our introduction of the concept of a disguised short burst leads to a 
third new important conclusion: the possibility that binary mergers can also lead 
to canonical long duration GRBs. It is clear, in fact, that if a binary system merges 
within a galaxy, then the average CBM density will be ^ 1 particle/cm^ and not 
the density typical of the halo (see Fig. [28l) . 

11. Ongoing progress in understanding GRBs, supernovae and 
neutron stars 

I would like now to address some ongoing research which promises to make signifi- 
cant progress on: 

(1) The analysis of the GRB-supernova connection, which is likely, among others, 
to lead to the first observations of a newly born neutron star, what I call a "neo 
neutron star" . 

(2) The discovery by the Fermi and Agile satellites of high energy emission from 
GRBs leading possibly to direct evidence of the ultra high energy baryonic 
component in GRBs. 

(3) The identification of the progenitor of the process of gravitational collapse; this 
has led to introduce the new concept of a "proto black hole." 

11.1. Progress in understanding the GRB-SN association 

I here recall a most interesting by product of the analysis of GRBs associated with 
supernovae. We examined GRB 980425 associated with SN iggSbw.ElIinHIIll in 
this system a peculiar long lasting emission in X-rays was identified and called 
URCA-1 in view of the fact that its first identification was presented at the Tenth 
Marcel Grossmann Meeting held in Rio de Janeiro .l^lSElSl \Yg also examined GRB 
030329, associated with SN 2003dh, clearly having a similar long lasting feature we 
named URCA-2.^^^ ^^^'^^^1 We have also studied GRB 31203, asso ciated with SN 
20031w, with a long lasting emission defined as URCA-3,E2MilIlIilsee Fig. W\\ and 
Tab. [5J I would like to review these results and present some general conclusions. 

11.1.1. GRB 980425 / SN 1998bw / URCA-1 

The best fit of the observational data of GRB 980425 we obtained with M.G. Bcrnar- 
^j^pMUnmil] ^^^^^ ^^ ^m ^ ^2 x lO^'s erg and B ^ 7.7 x IQ-^. This implies 
an initial e^e~ plasma with A'^e+e- — 3.6 x 10^^ and with an initial temperature 
T = 1.2 MeV. After the transparency point, the initial Lorentz gamma factor of the 
accelerated baryons is 7 = 124. The variability of the luminosity due to the inho- 
mogeneities of the CBlVp^is characterized by a density contrast Sn/n ^ 10^^ on a 
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length scale of A ^ 10^"* cm. We had determined the effective CBM parameters to 
be: (ricbm) = 2.5 x 10^^ particle/cm^ and (TZ) = 1.2 x 10~*, where TZ = Agff/Ayis 
is the ratio between the effective emitting area A^ff of the expanding shell and its 
entire visible area Ayis and takes into account both the effective porosity of the shell 
and the CBM filamentary structure! ^°^ * ^°^ l Details and fits of the data are presented 
in Refs. [2Tlti2T4l 

I then review the URCA-1 observations performed by BeppoSAX-NFI in the 
energy band 2-10 keVf2i2lby XMM-EPIC in the band 0.2-10 keVl^and by Chandra 
in the band 0.3-10 keV.'^SSlThe separations between the light curv es of G RB 980425 
in the 2-700 keV energy band, of SN 1998bw in the optical bandp2il22ll and of the 
above mentioned URCA-1 observations are given in Fig. |4TK. 



11.1.2. GRB 030329 / SN 2003dh / URCA-2 

For GRB 03032^213121612111 ^e obtained with M.G. Bernardini a total energy E^, 
2.12 X 10^^ erg and a baryon loading S = 4.8 x 10"'^. This implies an initial e+e" 
plasma with N^+^- — 1.1 x 10^^ and with an initial temperature T = 2.1 MeV. 
After the transparency point, the initial Lorentz gamma factor of the accelerated 
baryons is 7 = 206. The effective CBM parameters are {ricbm) = 2.0 particle/ cttt.^ 
and {TZ) = 2.8 x 10^^, with a density contrast 5n/n ^ 10 on a length scale of 
A ~ 10^^ cm. The resulting fit of the observations, both of the prompt phase and of 
the afterglow, were given in Ref. l216l217l We compare in Fig.HTB the light curves of 
GRB 030329 in the 2-400 keV energy band, of SN 2003dh in the op tical ba nd ^^^ l ^^^ l 
and of URCA-2 observed by XMM-EPIC in 2-10 keV energy band.I^SMUl 



^tot 



11.1.3. GRB 031203 / SN 2003lw / URGA-3 

The analysis of GRB 031203 with M.G. BernardinPsl2iMIll leads to a total energy 
El^ — 1.85 X lO^'' erg and to a baryon loading i? = 7.4 x 10^^. This imphes 
an initial e^e^ plasma with N^+^- = 3.0 x 10^^ and with an initial temperature 
r = 1.5 MeV. After the transparency point, the initial Lorentz gamma factor of 
the accelerated baryons is 7 = 132. The effective CBM parameters are {ricbm) = 
1.6 X 10^^ particle/cm'^ and (7?.) = 3.7 x 10^^, with a density contrast Sn/n ~ 10 
on a length scale of A ^ 10^^ cm. Particularly important in this source was the 
analysis of the luminosity in selected energy bands as well as of the instantaneous 
and time integrated spectra (see Figs. [40H2T|) . We recall here that such instantaneous 
spectra present a very clear hard-to-soft behavior, and that the corresponding time 
integrated spectrum is in very good agreement with the observed one. In Fig. I41C 
we compare the light curves of GRB 031203 in the 2-200 keV energy band, of SN 
20031w in the optical ban(pil222] ^^^^ ^j URCA-3 observed by XMM-EPIC in the 
0.2-10 keV energy bancp2^ and by Chandra in the 2-10 keV energy band.'22Sl 



July 6, 2011 0:9 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



63 



4.0x10* 



3.0x10' 



UJ 2.0x10* 



.::; 1.0x10* 



0.0x10" 



P-GRB 
position 



11-1+1 IlLlVj-Jftl T, 



'?f^^ 




2.5x10" 



2.0x10" 



1.5x10" 



1.0x10" 



5.0x10"° g 



0.0x10" 



-20 -10 10 20 

Detector arrival time (t^) (s) 



30 



40 



50 



-5.0x10" 



Fig. 40. Theoretically simulated light curve of the GRB 031203 prompt emission in the 20-200 
keV energy band (solid line) is compared with the observed data from Ref . 11151 The vertical bold 
line indicates the time position of P-GRB. See details in Ref. 11061 



11.1.4. Discussion 



In Table [2] (see details in Ref. I213P I have summarized the representative param- 
eters for the above four GRB-SN systems, including the very large kinetic energy 
observed in all SNe. 12221 Some general conclusions on these weak GRBs at low red- 
shift associated with SN Ib/c can be established on the grounds of our analysis: 

1) These results extend the applicability of our fireshell model to this low-energy 
GRB class at small cosmological redshift, which now spans a range of energy of six 
orders of magnitude from lO'^® to lO^"* erg! 93 | i°6 | i68 | 2ii | 2i6 | 2i7 |p^g^^^(,^|yg of ^j^jg ^.j^^gg 

is the very high value of the baryon loading which in one case (GRB 060218/^23 
is very close to the maximum limit compatible with the dynamical stability of 
the adiabatic optically thick acceleration phase of the GRBs.'^ Correspondingly, 
the maximum Lorentz gamma factors are systematically smaller than the ones of 
the more energetic GRBs at large cosmological distances. This in turn implies the 
smoothness of the observed light curves in the so-called "prompt phase." 

2) These four GRB sources present a large variability in their total energy: a factor 
10^ between GRB 980425 and GRB 030329. Remarkably, the SN emission both in 
their very high kinetic energy and in their bolometric energy appears to be almost 
constant, respectively 10^^ erg and 10*^ erg. 

3) The URCAs present a remarkably steady behavior around a "standard luminos- 
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Fig. 41. Theoretically computed light curves of GRB 980425 in the 2-700 keV band (A) P" I ^^^ I ^^ 
of GRB 03032 9 in the 2-400 keV band (B>2iMiL and of GRB 031203 in the 2-200 keV band 
^(2;'| 106 | 213 | 218 | ^j.g plotted, together with the URCA observational data and qualitative representa- 
tive curves for their emission, fitted with a power law followed by a n expon entially decaying part. 
The luminosity of the SNc in the 3000-24000 A range is also plottedBSHllllSee details in Ref. l2l3l 



ity" and a typical temporal evolution. We have considered the possibility that the 
URCAs are related to the SN event: either to dissipative processes in the SN ejecta, 
or to the formation of a neutron star in the SN explosion.'^i^ As an example, we 
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Table 2. 



GRB 


980425 


030329 


031203 


060218-^' 


(erg) 
B 


1.2 X 10**** 


2.1 X 10^2 


1.8 X 10^" 


1.8 X 10^'' 


7.7 X 10-3 


4.8 X 10-3 


7.4 X 10-3 


1.0 X 10-2 


70 


124 


206 


133 


99 


Tpholom 

(erg)'' 
(erg)b 


2.3 X 10-*'-* 


1.8 X lO^'S 


3.1 X 10*3 


9.2 X 10*** 
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Note: a) see Ref. 12271 b) Mazzali, P., private communication at 
MGll meeting in Berlin, July 2006; c) evaluated fitting the UR- 
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have given an estimate of the thermal emission from the neutron star surface and 
of the corresponding radius for URCA-1, URCA-2 and URCA-3 (see Table!!]). The 
different spectral properties of the GRBs and the URCAs have been pointed out.l^i^l 
It will certainly be interesting to compare and contrast the spectra of all URCAs 
in order to show the expected analogies among them. Observations of the URCA 
sources on time scales of 0.1-10 seconds would be highly desirable in presence of a 
strong enough X-ray signal. 

The investigation of the thermal evolution of neutron stars is a powerful tool 
to probe the inner composition of these objects. The cooling of neutron stars has 
been investigated by many authors, where many different microscopic models were 
assumed.'222H240| ]y[ost of the research on the thermal evolution of compact stars 
focus on objects with ages greater than 10-100 years, which is comprehensible if one 
considers that the thermal data currently available to us is for pulsars with estimated 
ages of or greater than 330 years. ^'^^^^'^^ The identification of the URCAs with young 
neutron stars, which we have recently called "neo neutron stars" (Negreiros et al., 
in preparation), will allow the study of the thermal evolution of young neutron 
stars in the little explored time window that spans from ages greater than 1 minute 
(just after the proto-neutron star regime^Sl) to ages < 10-100 years, when the 
neutron star becomes isothermal (see Ref. 12421 for details) . Relevant also are the 
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observations of the isolated Type Ic Supernova SN 1994p43land SN 2002aiPSI which 
show late emissions similar to the ones observed in URCA-1, URCA-2, and URCA- 
3. In our recent work with R. Negreiros we propose a revision of the boundary 
conditions usually employed in the thermal cooling theory of neutron stars, in order 
to match the proper conditions of the atmosphere at young ages. We also discuss 
the importance of the thermal processes taking place in the crust, which also have 
important effects on the initial stages of thermal evolution. 

There are three important ingredients that govern the thermal evolution of a 
compact star, these are: 1) the microscopic input, which accounts for the neutrino 
emissivities, specific heat and thermal conductivity; 2) the macroscopic structure 
of the star, namely its mass, radius, pressure profile, crust size, etc.; and 3) the 
boundary condition at the surface of the star, which provides a relationship be- 
tween the mantle temperature and that of the atmosphere, the latter being what 
we ultimately observe. These ingredients have been extensively studied, and a com- 
prehensive review can be found in Ref. 12351 As discussed in Ref. I242[ during the 
initial stages of thermal evolution (ages < 10-100 years), the core and the crust 
of the neutron star are thermally decoupled. This is due to the fact that the high 
density core is emitting neutrinos at a much higher rate than the crust, which causes 
it to cool down more quickly. This effectively means that initially the neutron star 
is cooling "inside out," with the core colder than the outer layers. This scenario 
is schematically depicted in Fig. 25] Details may be found in Negreiros et al., in 
preparation. 

All this suggests the exciting possibility that we are for the first time observing 
nascent, hot neutron stars. This possibility warrants further study in order to obtain 
a more concrete picture of the thermal evolution of neo-neutron stars. We have also 
encouraged dedicated observations of isolated SN in view of the similarities between 
URCA-l-URCA-3 and the Type Ic Supernova SN 199412^3 and SN 2002ap.[2iil 

11.2. The high-energy emission in the fireshell scenario 

The launch of the AGILE and Fermi satellites opened a new high energy window for 
GRB science. Thanks to the Large Area Telescope (LATp^ on-board the Fermi 
spacecraft, it is possible to detect very high energy photons (> 100 MeV) from 
GRBs. The first results obtained from Fermi is that not all GRBs have this high 
energy component, but just a fraction of them: currently LAT has detected energetic 
photons from ^ 20 GRBs among the whole set of ~ 480 GRBs detected by the Fermi 
Gamma-Ray Burst Monitor (GBM). '233 However, this last evidence has already been 
noticed by the EGRET detector on board the Compton Gamma-Ray Observer 
satellite, which detected high energy photons only from a few GRBs. 

The other very interesting discovery of the Fermi LAT was that this high energy 
component is delayed compared to the emission observed by the GBM detector. 
Moreover when the GBM signal fades in the canonical way, the high energy one 
is prolonged and can last for hundreds of seconds. The case of GRB 080916cf^ is 
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Fig. 42. Schematic representation of the cooling of a young neutron star. Due to stronger neutrino 
emissivities, the core of the star cools down more quickly than the crust, causing the star to cool 
inside out. Darker and lighter areas represent higher and lower temperatures respectively. Details 
in Negreiros et al., in preparation. 



remarkable (see Fig. 143)) : it is currently the most energetic GRB to date, with an 
estimated isotropic energy of Eiso = 8.8 x 10^^ erg. A significant contribution to 
this value was provided by the high-energy gamma-rays (> fOO MeV) emitted by 
this GRB. Moreover this high energy emission started ~ 3.6 s after the GBM trigger 
and persisted up to 1400 s after the GBM trigger. The spectral energy distribution 
of these high energy photons is best represented by a power-law function whose 
photon index was found to be the same as the harder component of the GBM Band 
spectrum. This observed continuity does not appear to be an universal property of 
all sources presenting a high energy emission. 

The main contribution of the fireshell model applied to this family of sources is: 

1) That the initial phase of GRBs lacking the high energy component represents 
the P-GRB. 

2) The corresponding determination of the Lorentz gamma factor of the source, 
and the strong correlation between the existence of the high energy component and 
very large values of this factor (7 > 10'^). 

3) The possible connection of the high energy component to the ultrarelativistic 
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Fig. 43. Fermi GBM and LAT light curves of GRB 080916C. Details in Ref.[247] Note the P-GRB 
emission starting before the LAT data. 



baryons interacting, via proton-proton collision, with a dense CBM following the 
classical work of Fermi (see e.g. Ref. I248p . 

I briefly recall some preliminary results on the parameters of GRB 080916c 
within the fireshell model. The best fit correspond to a value Edya = 8.8 x 10^^ erg. 
From the observed flux, we estimated the energy emitted in the P-GRB, Ep-grb = 
2.8 X 10^"^ erg, which, in view of Fig. [TH gives a value of -B = 3.3 x 10^^. The 
corresponding value of the average CBM density is {n) — 2.2 x 10"'^ particles/cm~^. 
From these values, we infer a very high value for the Lorentz gamma factor at the 
transparency: 70 = 3.17 x 10'^. 

The identification of the first pulse in GRB 080916c with the P-GRB is further 
supported by its duration: it is comparable with what is expected for a P-GRB, e.g. 
on the order of few seconds.^ Moreover, a time-resolved spectral analysis of the P- 
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Fig. 44. Spectral fits of the Fermi GBM and BGO count spectra of the P-GRBs of GRB 080916C 
in the range (8 keV— 5 MeV). We consider in these cases a blackbody spectral model with an 
extra power-law component, a classical band function and a modified black body spectral model 
respectively. 



GRB appears to be consistent either with a classical band model, represented by two 
broken power-laws smoothly connected at a given energy Eq or with a black-body 
model with an extra power-law component (see Fig. I44p . 

All the above represents a first step toward the understanding of these GRBs 
within the fireshell model and opens further studies on the origin of the high energy 
component by refinements of the energy release process in the comoving frame of 
the baryon component interacting with the CBM. 



11.3. Possible multiple components in a GRB: the case of GRB 
090618 

I have emphasized how the study of GRBs encompasses an unprecedented number of 
fundamental physics topics, each one probing untested regimes. The physics of grav- 
itational collapse leading to the formation of a black hole needs a new self-consistent 
treatment in the ultra-relativistic regimes of strong, weak, electromagnetic and 
gravitational interactions as pointed out in Sec. 14.21 The outcome of the collapse 
leads to the formation of an optically thick e'^e^-baryon plasma, self-accelerating 
to Lorentz gamma factors in the range 200 < 7 < 3000, never reached before in 
macroscopic objects. These novel ultrarelativistic plasma regimes have been sum- 
marized in Sec. 14. II The transparency condition of the e"'"e~-baryon plasma, already 
observed in the cosmological background radiation, gains in GRBs the additional 
component of ~ iq56-57 ^j^ryons in the TeV region (see Sec. [3]). The collision of 



15-16 



cm, 



these baryons with the CBM clouds, characterized by dimensions of 10 
leads to high energy collision regimes never reached before in dimensions and ener- 
getics (see Sec. 111. 2() . Before concluding this review, I would like to outline how the 
progress in the last decades of optical, X and gamma-ray observatories has made 
possible the observation of GRB 090618 offering an unprecedented possibility for 
testing the theoretical models of the above objects with crucial observations and 
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Fig. 45. Evolution of kT (left) and power-law index s (right). The grey curves are the detector 
count light-curves, with arbitrary normalization. The time is measured from the trigger, except 
in the plot of kT for GRB950624, where time is measured from the beginning of that pulse. See 
details in Ref. \7ET\ 



linking for the first time the observation of GRBs to their immediate progenitor. 

GRB 090618 is indeed one of the closest (z — 0.54) and most energetic {Eiso = 
2.70 X lO^'^ erg) GRBs, observed under ideal conditions by the satellites Fermi, 
Swift, Konus-WIND, AGILE, RT-2 and Suzaku, as well as from on-ground optical 
observatories. It therefore represents an ideal case to test our fireshell model. We 
have analyzed the emission from this GRB with special attention to the thermal and 
power-law components highlighted by F. Ryde and collaborators. It is interesting in 
fact that independently of the development of new missions, the BATSE data have 
continued to attract full scientific interest, even so many years after the end of the 
mission in 2000. Based on the BATSE data, important inferences for the spectra 
of the early emission of the GRB have been made in Refs. 12491 and [2501 They have 
convincingly demonstrated that the spectral feature composed of a black body and 
a power-law behavior together plays an important role in selected episodes in the 
early part of the GRB emission. They have also shown in some cases a power-law 
variation of the thermal component as a function of time, following a t^" behavior 
(see Fig. [45]) . 

We have proceeded to the determination of the fundamental parameters of GRB 
090618 in the fireshell model, including the identification of the P-GRB, the dyado- 
sphere energy, the baryon loading, the density and porosity of the CBM. We have 
clearly shown the existence of two different components in GRB 090618: Episode 
1 and Episode 2 (see Fig. [46]) . We have confirmed that Episode 1, lasting 50 s, 
cannot be considered either a GRB or a part of a GRB .1222] Much like the cases 
considered by F. Ryde, it shows a very clear thermal component with temperature 
evolving following a precise power-law behavior between kT — 54 keV and kT = 12 
keV (Izzo et al., in preparation). We have proposed that the radius is increasing 
between ^ 12000 km and 70000 km, with an estimate mass of ~ 15M0. The second 
component is a canonical long GRB with a Lorentz gamma factor at transparency 
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Fig. 46. The two episode nature of GRB 090618. Details in Ref. [253l 



of 7 = 490, a temperature at transparency of 25.48 keV and with characteristic 
size of the clouds, generating the observed luminosity variations, of Rd ~ 10^^ cm 
(see Fig. |47]). Episode 1, in our interpretation, is related to the progenitor of the 
collapsing bare core leading to the black hole formation: what we have defined a 
"proto-black hole" of ~ 15Mq. For the first time we may be witnessing the pro- 
cess of formation of the black hole from the phases just preceding the gravitational 
collapse all the way up to the GRB emission. 

Of particular interest in this respect are the first two-dimensional hydrodynam- 
ical simulations of the progenitor evolution of a 23Mq star close to core-collapse, 
leading to a naked core, as shown in the recent work of Arnett and Meakin.'^SSl 
In that work, pronounced asymmetries and strong dynamical interactions between 
burning shells are seen: the dynamical behavior proceeds to large amplitudes, en- 
larging deviations from the spherical symmetry in the burning shells. It is of clear 
interest to find a possible connection between the proto-black hole concept that we 
have introduced and the Arnett and Meakin results: to compare the radius, the 
temperature and the dynamics of the core that we have found with the naked core 
obtained by Arnett and Meakin from the thermonuclear evolution of the progeni- 
tor star. There is also the interesting possibility that the CBM clouds observed in 
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Fig. 47. Simulated light curve of the extended afterglow of GRB 090618. See details in Ref. [2521 



GRBs, as well as the baryon loading captured by the e^e^-plasma dynamics, are 
related to the vigorous dynamics in the violent activity of matter ejected in the 
evolution of the original massive star prior to the formation of the naked core. 



12. Conclusions 

I have summarized the vast scenario made possible by the unprecedented techno- 
logical possibilities of observing SNe, GRBs, GRBs associated with SNe, all the way 
to neo-neutron stars and proto-black holes. The corresponding theoretical work is 
rapidly expanding over new regimes of ultrarelativistic field theories and encom- 
passing all fundamental interactions. A novel unified approach to nuclear physics, 
neutron stars, white dwarfs and gravitationally collapsed objects is rapidly emerg- 
ing. All this offers unprecedented possibilities to enlarge our knowledge of funda- 
mental physics and relativistic astrophysics, as well as to reach an understanding 
of yet unexplained phenomena of fundamental physics. 
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